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Methodology of including QCD NLO corrections in the quark-gluon Monte Carlo shower is outlined. The 
work concentrates on two issues: (i) constructing leading order (LO) parton shower Monte Carlo from scratch, 
such that it rigorously extends collinear factorization into the exclusive (fully unintegrated) one which we call the 
Monte Carlo factorization scheme; (ii) introducing next-to-leading-order (NLO) corrections to the hard process 
in this new environment. The presented solution is designed to be extended to the full NLO level Monte Carlo, 
including NLO corrections not only in the hard process but in the whole shower. The issue of the difference 
between the factorization scheme implemented in the Monte Carlo (MC) solution and the standard MS scheme 
is addressed. The principal MC implementation is designed for the electroweak boson production process at the 
LHC, but in order to discuss universality - process independence, the deep inelastic lepton-hadron scaterring is 
also brought into the MC framework. 
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I. INTRODUCTION 

The excellent performance and fast experimental data accu- 
mulation of the Large Hadron Collider (LHC) at CERN makes 
the precise evaluation of the strong interactions effects within 
perturbative Quantum Chromodynamics (QCD) |fl}{3] a more 
and more important task. The principal role of QCD in hadron 
colliders data analyses (LHC and Tevatron) is to provide pre- 
cise predictions for distributions and luminosities of quarks 
and gluons accompanying production of heavy particles. 

One of the most important theoretical tools of perturbative 
QCD (pQCD) are factorization theorems |@]-0, which refor- 
mulate any scattering process in QCD in terms of the on- 
shell hard process part convoluted in the lightcone variable 
with the ladder parts, provided a single large scale Q 2 is in- 
volved (short distance interaction). The hard process is usu- 
ally treated at a fixed perturbative order and the ladder parts 
are resummed to infinite order, for each coloured energetic in- 
going/outgoing parton. The initial state ladders give rise to 
PDFs, the inclusive parton distribution functions. 
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The initial state ladder, instead of being source of the in- 
clusive PDF, can be also modelled using Monte Carlo sim- 
ulation (including hadronization), as initiated in refs. |7][8]. 
Such an implementation of the QCD ladder is referred to as 
the parton shower Monte Carlo (PSMC) program. This kind 
of programs play an enormous practical role in all collider ex- 
periments. In today PSMCs the initial ladders are restricted 
to the leading order (LO). With growing requirements on the 
quality and precision of the pQCD predictions for the LHC 
experiments it becomes urgent to upgrade PSMC to the same 
next-to-leading-order (NLO) level, which was reached for the 
inclusive PDFs two decades ago. This is not easy, mainly be- 
cause factorization theorems of QCD |4-6| were never meant 
for the MC implementation. They are well suited for the sim- 
pler case of the hard process upgraded with the finite order 
calculations convoluted with the collinear inclusive PDFs. 

There has been, however, a significant progress in imple- 
menting pQCD in the framework of PSMC, which started with 
the work of MC@NLO team |9l[T0]], and followed by develop- 
ment of POWHEG method lfTTl[T2"l . In these works hard pro- 
cess in PSMC is upgraded to the NLO, while the ladder part 
stays at the LO level, essentially older solutions and software 
for the LO PSMC are not modified. This, of course, saves 
a lot of work but because of that the methodology of com- 
bining the initial ladder parts and the NLO-corrected matrix 
element (ME) for the hard process is quite complicated. The 
solution for this problem is to redesign the basic LO PSMC. 
This would be too big an investment, if simplification of the 
NLO corrections to the hard process ME were the only aim. 
However, this effort is mandatory if we are aiming to upgrade 
also the ladder parts of PSMC to the complete NLO level. 

In this paper we outline a redesigned LO parton shower MC 
and simultaneously present methodology of including QCD 
NLO corrections to the hard process which takes advantage of 
it. Hence we shall concentrate on two issues: (i) constructing 
once again the LO parton shower Monte Carlo from scratch, 
such that it is based firmly on the rigorous extension of the 
collinear factorization theorems, which, contrary to the orig- 
inal collinear factorization, is fully exclusive (unintegrated), 
(ii) introducing the NLO corrections to the hard process in this 
new environment. It is natural to expect that the issues of the 
difference between the factorization scheme implemented in 
the Monte Carlo (MC) solution and the standard MS will have 
to be addressed. The important point of universality (process 
independence) will be also discussed extensively. Although 
the principal aim will be a new MC implementation with the 
LO ladder (upgradable to NLO) and the NLO hard process 
for the production of the electroweak (EW) boson in quark- 
antiquark annihilation, in order to address the issue of uni- 
versality, and for other practical reasons, the deep inelastic 
lepton-hadron process will be also brought into the consider- 
ation. The next step in the project, that is the upgrade of the 
ladder part to the NLO level, will be treated in a separate pub- 
lication 1131 . although the general method was already out- 
lined in refs. Ifl4l - [l6l . Many technical details needed for the 
NLO ladder are provided in ref. ifTTl and auxiliary discussions 
on the soft limit and the choice of the factorization scale in the 
MC can be found in refs. 1 18-201. The first numerical tests of 



the discussed method, of including the NLO hard process cor- 
rections, are presented in ref. GTl . 

Let us note that the ongoing effort undertaken in refs. IT221 
23 1 is in some aspects similar to the present work, in particular 
the parton shower is also redesigned at the NLO level and the 
departure from the standard MS factorization scheme is also 
advocated. These works are extending/exploiting techniques 
of refs. El E21. 

The remaining part of the paper is organized as follows. In 
sec. |ll]we discuss collinear factorization in the form suitable 
for the MC implementation. Section III covers construction 
of the LO MC for EW boson production and the new method 
of introducing the NLO corrections to the hard part of this 
process. In sec. IV we present a similar MC solution for the 
DIS process, with the new LO MC modelling of the initial and 
final state ladders and the NLO corrections to the hard process. 
The issue of the universality as well as factorization-scheme 
dependence are addressed in various steps of this presentation, 
with the final discussion in sec.|V] where we also summarize 
and give outlook of further work. 



II. GENERALITIES - COLLINEAR FACTORIZATION 




FIG. 1 . The EGMPR |4) factorization for the EW boson production. 
Example 2PI kernels for a quark in the forward hemisphere Kof , 
an antiquark in the backward hemisphere A'yg and the hard process 
part Co are delimited by ellipses. The lower figure highlights the use 
of cut diagram notation by indicating the amplitude in blue and the 
complex conjugated amplitude in red. 



The precise definition of the LO approximation within fac- 
torization of collinear singularities, and the resulting distribu- 
tions implemented in the Monte Carlo, is a necessary prelude 
of defining complete NLO distributions, both in the ladder and 
in the hard process. This is why in the following we shall de- 
fine a new LO MC, "anchoring" it in the collinear factoriza- 
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tion theorems |4-6| as firmly as we can. 



The production process of electroweak (EW) bosons W, Z, 
y, in the hadron-hadron scattering, or of other colour-neutral 
vector particle, will be the object of the interest in the fol- 
lowing. We shall refer to it as Drell-Yan (DY) process for 
short. We are going to describe a Monte Carlo algorithm in 
which two initial state parton ladders will be modelled up to 
the LO and the hard process to the NLO level, both the hard 
and ladder parts in a completely exclusive way. In the con- 
struction of the MC we will keep track of the precise rela- 
tion to the QCD factorization theorems, having in mind that 
the ladder part (parton shower) will be upgraded to the NLO 
level in the next step. Diagrammatically, we shall temporarily 
limit ourselves to the Cp part in the ladder, that is to gluon 
bremsstrahlung. Because of copious soft gluon production, 
this is the most difficult part in the MC construction of the 
LO (later on of the NLO) MC implementation for the lad- 
der. Adding more diagrams (quark-gluon transitions, singlet 
diagrams) will be discussed briefly, but will be treated in a 
separate publication. 



Let us start from the "raw collinear factorization" formula 
of ref. [4 1 (in the axial gauge) illustrated in Fig.[T]in a standard 
way (cut diagrams). Following closely notation of ref. 0], the 
standard Feynman amplitude for the heavy boson production 
process is ^j F fJ\PF ,p\B,<7i,<72,r), where two incoming par- 
tons of the type jp and jg have spin indices rp and rg, heavy 
boson decay lepton momenta are qi and any number of the 
emitted on-shell gluons and quarks are collectively denoted 
as r. Denoting also collectively (pFtpB,q\,<l2) as (p,q), the 
partly integrated cut-diagram is defined as 



r J 



xAy;? B { P , q ,r)A%i(p, q ,r) 



(i) 



It is related directly to differential cross section daj F j B (p,q) — 
jWx^r F r B ^j F j B rFrB (p,q) Z F Z B , where Z FM are wave renor- 
malization factors. The essence of the "raw" factorization the- 
orem in ref. |4| is that all collinear singularities are located in 
the ladders with multiple 2-particle irreducible (2PI) kernels 
K r J, r J(k,p). Suppressing for the moment ladder B and ne- 



glecting subscript F the above statement reads 

MJ s (p,q)=CoJ(p,q) + Y,fl [d 4 Pm 

X £ C^»(p n ,q)K Q ^ ,S "- l (Pn,Pn-l) 
X L K 0k-^jn-2 (Pn-l,Pn-2)- 

r n-l 



Using the compact matrix notation of refs. [4, 26], the above 
expression in case of two ladders reads 



O = C 



1 



1 



1 — -K p 1 — -K 0B 



= e Co (-wow 



n \ ,«2— 



(2) 

see also the upper part of Fig. [T]for an equivalent graphical 
representation. In the lower part of Fig. [T] the above for- 
mula is illustrated diagrammatically using the lowest order 
bremsstrahlung matrix element where we explicitly indicate 
the 2PI kernels, with the red part of the diagram representing 
the conjugate part A* of eq. ([TJ. Note that in the above expres- 
sions phase space of the emitted on-shell partons (cut lines) is 
integrated over and treated inclusively. In the following dis- 
cussion it will be explicit and implemented in the Monte Carlo 
parton shower. 

According to ref. |4| all collinear singularities in o of eq. |2]) 
are coming from (dressed) propagators between kernels Kqf 
{Kqb) along the ladders. The 2PI kernels for the initial quark 
ladder, Kqp, and the antiquark ladder, Kqb, are expanded to the 
infinite order, see ref. (4|. In the following practical example 
we shall truncate them to the (lowest) first order (LO) or to the 
second order (NLO) K^p = K^p +Kqp, taking into account the 
following ~ Cp diagrams: 



K, 



(i) 

OF 



K, 



OF 



K, 



[2] _< 



OF 




A similar expansion up to the 1st order (NLO) is done for 

the hard process part Cq 1 ^ = C °' +Cq . For simplicity, we 
are omitting the initial quark and antiquark distributions in 
the beam hadrons and the flux factor is included in Co- The 
dot "•" in the product A ■ B means full phase space integration 
/ d 4 q over the lines joining two subgraphs in one ladder. 

The next step in the classic works of refs. [4, 26] is the in- 
troduction of the projection operator P. Its role is to decouple 
kinematically not only Co and the ladder parts, but the consec- 
utive kernels Kq along the ladders, such that the integration 
over lightcone variables and collinear logs become manifest 
and ready for analytical calculations. Formally eq. (|2]) gets 
transformed (at the infinite order) into 
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11°° 

o = c®r F ®r B = c ——— ———= V c((g>K F )"i(®K B ) 

1-®K F 1-®K B nin2=Q 
1 1 1 

C — Cq — — — — — , Kp — — — j , Kg 



1 



(3) 



R=(-1-®P), 



1 - ®TPK 0F 



l 

i-hk of 



1 - ®PK 0B 



l 



l-BJSbB 



where A (8 B means convolution in the lightcone variable 
/ ifzicfeSO* — ziZ2M(zi)B(z2)> while integration over trans- 
verse momenta is traded into \ poles of dimensional regu- 
larization, extracted (upon integration) by P, see ref. [26| for 
details. The projection operator P used in is slightly dif- 
ferent. However, both approaches are incompatible with any 
MC implementation, as can be seen from the explicit expan- 
sion up to NLCQ 



C = C?>(1 



OF 



OB) 



OF ' 



V 0F 



[1] 



of ' a of) 



(4) 



,P<) + (P<) 



Why? As we see, the above is a mixture of the original phase 
space integrals like (Kq F -Kq F ) and of partly integrated inte- 
grals like in (WKo F ) ® (WKq F ). Even if we managed some- 
how to undo the transverse momentum integrations implicit 
in P operator, we would still face huge (double logarith- 
mic) oversubtraction in ~PKq F ((-1 — ®¥)Kq F ) compensated 



by (PK [ q F ) ® (P^), which would be deadly for any MC 
implementation. For an explicit demonstration of this prob- 
lem see also the toy model considerations in ref. |27||^]or LO 
analysis to infinite order in ref. (16). 



qq-*-Z+nG 




p, ^fvnpi-y/^m 
p, fknj?kry( 



FIG. 2. The EGMPR factorization for EW boson production. The 
LO ladder with the projection operators P' inserted between the LO 
kernels. 

The solution of the above over-subtraction problem is 
well known and already employed in the existing LO parton 



1 Omitting for simplicity quark wave function renormalization. 

2 In this work another example of the operator P is presented and also the 
order in factorizing collinear singularities in eqs. J3} is reversed - it starts 
from the hard process. Nevertheless, it features the same oversubtraction 
problems that inhibits MC implementation. 



shower MC since long Q H) - in short, one has to introduce 
the time-ordered (TO.) exponent, see also ref. [28|. Beyond 
the LO, a collinear factorization formula with TO. exponent 
was outlined in ref. ifTBI and we shall adopt it here. We are go- 
ing to use it for the LO MC ladders combined with the NLO 
hard process ME0 According to ref. [16) eq. ^ is replaced 
by: 



n\ — 1 «2 — 1 



cS°){(p<)r} r . a {(p'^r} r . a , 



(5) 



where K^J is the lowest order (LO) 2PI kernel, the same as in 
eq. ([3]), but at the NLO and beyond it is different, see the defi- 
nition in ref. IfTBI . The above LO process is depicted in Fig. [2] 
As shown in ref. [16], the complete and rigorous definition of 
the new projection operator P' is not simple. For the present 
purpose of the two LO ladders and the NLO hard process in 
the EW boson production we shall define it step by step, start- 
ing from the simple cases of zero, one and two gluon ME, 
n\ -\-n% = 0, 1, 2, rather than defining it immediately in the full 
form. 

Let us denote the Born level differential cross section for 
the EW vector boson production and decay process, qq — > 11, 
as follows: 



d£l 



It may also include non-photonic EW radiative corrections. 
The above differential cross section is so well known that we 
may avoid defining its details explicitly. From (two) Feynman 
diagrams we obtain the following (exact) LO single-gluon 
emission differential distribution, 



dot = 



C F a s dadfi dq 

do B (s,Q F ) (1-p) 2 
dQ. 2 



do B ($,Q B ) (l-a) 2 
dQ 2 



dQ., 



(6) 



where the Sudakov variables a and p are defined in Ap- 
pendix [A] This elegant formula is valid for any on/off-shell 
vector particle production, B = y,W,Z. The polar angles Q FB 
are defined |29| with respect to — po B and po F respectively in 
the rest frame of the B boson (the rest frame of po F + po B — k). 



3 Considerations concerning the NLO MC ladder can be found in refs. 1151 - 
[17). 
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Two collinear limits p — !• 0, 1 — a = z = const or a — >• 0, 
1 — p = z = const in eq. ^ are manifest. For instance, in the 
1st case we have: 



d<3] = 



C F a s dzdfi dy 1+z 2 da B (zs,Q F ) 



(l-z)P 2ti 2 



da 



d£i, (7) 



that is the LO kernel P qq (z) = ^j— 4yz^j shows up as ex- 
pected. Introducing the P' projector in this context may look 
like an overkill, but it will be instructive to explain how it 
works in this simple case before going to not so obvious case 
of multiple use of P' in the following. 

The necessary ingredient is a spin projection operator P vpi „, 
which we define a little bit more rigorously as compared to 
the Curci-Furmanski-Petronzio work (CFP) |26|. Our P sm „ 



(8) 



acts definitely before the phase space integratioij^] 



HP spin K = H4 t 



An ■ qt 



where qi is the on-shell momentum entering the H part, such 
that it conserves the longitudinal (lightcone) component n ■ 
q~j = n- qi, and for the axial gauge vector n we may take n — 
Pob or any other light-like vector, the same for all rungs in a 
given ladder. The same n is defining transverse polarizations 
of gluons in the axial gauge and enters a definition of lightcone 
variables Xj = (nqi) / (npop). 

In the present case of a single gluon emission, the action of 
the P' projector, inserted in the squared, spin summed Feyn- 
man diagrarrj^Jbetween the gluon emission vertex and the qqB 
vertex (B is the EW vector boson) can be summarized as fol- 
lows: 

(a) Apply P sp i„ to decouple a spinor y-trace into two parts, the 
hard ME part and the ladder part. 

(b) Apply the explicit upper limit of the phase space for an 
emitted gluon in the transverse momentum, for example 
a <M, with a 2 = s$/ a. 

(c) Take the expression for the hard part MB^jin the collinear 
limit, a — > (p —> 0) keeping a = 1 — z = const (or a + P = 
1 — z = const) and extrapolate it all over the phase space. 

(d) Keep unchanged the phase space integration element and 
its limits. 

Point (c) of the above recipe is the most important and requires 
more discussion. Finding out the limiting collinear expression 



is trivial, see eq. |7). What is non-trivial is the off-collinear 
extrapolation (OCEX) of this formula, out of p = point to 
all non-collinear phase space. The simplest recipe: go back 
along z = 1 — OL line and use eq. (|7]i 



daw = 



Oa s dadfi afcp 
n ap 2k 

l + (l-a) 2 da B ((l 



a)s,B f 



(9) 



da 



da, 



would be acceptable, provided the Born cross section is flat. 
In the presence of a narrow resonance in dogis), in the s = 
s(l — a — p) = sz variable, this would lead to a disastrous NLO 
correction dO\ — dOw, wildly varying over the phase space. 
This kind of OCEX follows a vertical dashed line in Fig.|3ja). 

However, there is a freedom in the off-collinear extrapola- 
tion away from p = point - we may do it also along the line 

x = 1— a — p = const: 



dOLO = 



C F a s dadfi c/cp 
7t ap 2n 

l + (l-a-P) 2 £/a B ((l-a-p>,e) 

2 d~a 



(10) 



da. 



In Fig.[3|a) this kind of OCEX goes along the curved dashed 
line there. The angle 8 has to be also defined within OCEX, 
in any reasonable way which coincides with the correct value 
Qf at the P = poinQ In a sense, the above is fully com- 
patible with the methodology of calculating the NLO correc- 
tions to the EW production process in ref. Il30l where a light- 
cone variable of the collinear factorization is also mapped into 
x = 1 — a — p. The essential difference is that with the help 
of OCEX we are replacing the traditional collinear PDF of 
ref. (30j with the exclusive distribution of eq. JT0j>. A number 
of consequences of this replacement will unfold gradually in 
the following. 

For easier generalization to two and more gluons, let us for- 
malize slightly the above as: 

^ M = PPa6M>a Pspin, 

where PT^ a takes collinear limit and implements off-collinear 
extrapolation on both sides of its location in the Feynman dia- 
grams, without affecting the phase-space integration element: 



4 In CFP spin projection is part of HP K = -f/^,] P E [ K, where P e ex- 
tracts the pole part and simultaneously sets q 2 — > in H. We have to be 
more specific about the choice of q in the substitution q — > q, q 2 = 0. 

5 Just one Feynman diagram for gluon emission from the initial quark. 

6 The one-gluon ladder can remain unchanged, but in the case of two gluons 



in Kq taking a limiting expression is also done in the ladder. 
7 For instance, as an angle between qj — qi and pof — poe in the rest frame 
of the EW boson. 
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► lnoc 




► lna 



(a) (b) 

FIG. 3. The Sudakov plane of single gluon emission in DY for angular ordering. 



cfpj,*$ = /^4fW P M JdT 2 (P-k;q u q 2 )C { ° ] (q u q 2 ,k) 
2CfOL, dadqdu l + (l-a) 2 cIgb(s(1 - a). 8) 



(11) 



r 



Note that z = 1 — OC = i/s is within the proper limits < 
z < 1 . The above can be formalized by means of introducing 
the rescaled 4-momentum W = X^ l k^, 



a 



a 



< 1, 



where a = a/A, and P = p/X. The dilatation transformation 
preserves angles, hence in the MC one may generate a and 
a = 1 — z according to 

Co^'m K QF = J ~jjX) ^ K 0F (*) P M 



/ 



d% 2 {P-k;q u q 2 ) C^\q\ 1 q 2 ,k) 



2C F a s dady dz 1+z 2 Ir . da B {sz,Q) a 

o a <M ^ dil q — o lz >0, 



a 2n 1 — z 2 



(12) 



then construct W and finally rescale it k = Xk. The 
"barred space" of W depicted in Fig. |3jb) is merely a re- 
parametrization of the true phase space in Fig. |3ja). For a 
single gluon the above parametrization of the phase space may 
look trivial, but for many gluons it will be useful. 

As already said, the main role of the P operator in refs. |4| 
l26l is to decouple kinematically the hard process and the lad- 
der. The above P' does it also, but more gently, protecting the 
4-momentum conservation. The kinematic decoupling is seen 
from the phase-space integration 



r (0)p/ 



M A 0F 



ln- 



M 
i — 

qo 



l 



1 , 2C F a s 
o dZ —W 



-z) 



which provides exactly the same result as the collinear factor- 
ization in these classic works where the 4-momentum conser- 
vation is broken. An additional cut-off a > qo (phase space 
boundary in the LO MC) was used in the above. 

One may also easily define, within the above scheme, a pro- 
totype of an universal inclusive collinear PDF: 



D(M,x) 



In 



OF 

M 2C F a s 



(14) 



qo 



2(1- 



thanks to P^ closing the phase space from the above by 
means of the factorization scale M. As we shall see later on, 
in the analogous construction for the deep inelastic electron- 
proton scattering, the role of W' M is to map the phase space 
of the ladder into an idealized phase space of kf, decoupled 
kinematically from the hard process, thus removing process 
dependencejjand gaining universality of the ladder part! 



A. NLO correction to hard process - one real gluon 

Having defined the sing le-gluon ladder parts cjp~P'K$ 
and Cq°'p'Xqj in the exclusive way, within the same exact 



<3g (sx) , (13) 8 ^ phage space shape (upper limits) usually depends on the process type. 
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phase space where the complete exact 1-gluon distributiorj^] 
of eq. d6}, 



is defined, it is straightforward to define the LO + NLO factor- 
ized hard process part in the exclusive (unintegrated) manner: 



-Q 



— L o 

-c [( 



i + (i-p')#( (1 



OF 



(1-P')^, (1) 



OB 



■cfi+c{'(i-p')«& , +cfi(i-p')^ 



■ck 



r [0]jAi) 

-«- OF 



r [0]^(l) 
"°0 OB 



OB 

-C n ir A 0f . — C n Jr A nD 



V 0B 



.c^p'^). 



(16) 



D [i] -r [1] +C l0] -K W - 

U — U +U Q "-OF 



[0] ^(i) 



OB 



The above difference of the exact and approximate ME at the 
(15) level of the integrand reads 



C 11 +C ( ? )1 (l-P')4F ) +C ( f 1 (l-P')^B 



/ 



— ^ dLl a 

ap 2ti ' 7C 



(i 

B 

(1-p) 2 da B 
2 d£l a 



(f,6f 



(i) 

(1 - a) 2 afa B 

~~2 



(f,e fl ) 



(17) 



l + (l-a-P) 2 rfc B fl i + (i- a _p)2 dCT 



dxC2r{x) Gb(sx) 



where the C 2r {x) function is calculated in the Appendix [b] 

Note that in the above integral the phase space is covered 
by the LO distribution completely, without any gap or over- 
lap, provided the factorization scale variable M in both hemi- 
spheres is adjusted conveniently. Similarly, the entire inte- 
grand of the NLO correction is also defined all over the 1- 
gluon phase space. 



B. LO two-gluon C^P'K^P'K^ ■ 
ladder MC 



a prelude for the LO 



The factor A. 2 from the phase space is compensated by a sim- 
ilar factor in the 1-gluon distribution, as it is for the 1-gluon 
case. No approximations nor projections are present yet. 



Insertion of the first P' requires examination of the 
collinear limit a\ — >• (pi — s- 0) keeping oti = const; in this 
limit we have also pi — s- 0, OL\ — > a and Xi—^l, The spin pro- 
jection operator of eq. (|8]l is also used. The collinear limit is 
well known 



Let us take the exact integrated distribution for the ladder 
diagram (no projections) with two gluons emission 



2 = C, 



(o) AV) 

'^OF -A 0F 



d 3 kid 3 k 2 , . , 

dx\ — n Trdiiyr — k\ 

2k\ 2k\ y 



PB(poF,poB,hM,qi,q2) s8 sx , P _ k 



■k 2 ;qi,q 2 ) 



(18) 



in which we have introduced explicitly a variable being the 
effective mass squared 



fl =sxi = (q\ + q 2 - 



-k 2 f = {P-k l f 



of the final state system after emitting the gluon k\ . 

Let us start with the same operation of parametrization of 
the phase space in terms of k\ = X\k\ as in the previous case 
of the single gluon, for a gluon at the end of the ladder: 



r (oMiMi) _ 

<-0 -^OF^OF ~~ 



, d 3 hd 3 k 2 ^ 



h = 



x 'KiPB(poF,poB,h>ih,k 2 ,qi,q 2 ) B^i-bj, 
s(l-xi) _ ai 



2Pk, 



oci+Pi 



(19) 



2 T 1 2C F a s l + (l-ai) 2 
hm a\ Pa(Afei,*2)—) = — j 



x 



2 (20) 
Pi(xipoF,PoB;k 2 ,qi,q 2 ), 



where pi is the already discussed distribution of the single- 
gluon emission p of eq. (|6j, in the reduced center-of-mass 
system of sx\ — (xiPof + Pqb) 2 , provided we rename k —> k 2 . 



The formula obtained above, valid originally at the 
collinear point a\ = 0, is now extrapolated to the off-collinear 
phase space using k\: 
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a 2 F - <- K 0F if K QF 



dx 



d?k\ 2 d i k2 



2k' 



q-A-! — -Q-di2{P-'kiki -k 2 ;qi,q 2 ) 



2C F a s 



2k 1 



3a2>ci[ 



l + (l-OCl 

2afo^ 



pi(xipoF,poB;k 2 ,quq 2 )8 Xi 



= l-OL| 



(21) 



dxi8 x . 



Cf(X s P(zi) da\ c/cpi 
7t 1 — zi a\ 2% 



Qa 2 >a i dlj,(P - r k\k\;qi 1 q 2 M) 9\{z\Pqf ,P0B;k 2 ,qi,q 2 ) , 



where the LO splitting kernel for the first emission of k\ 



4h) 



P(z) _ 1+z 2 
1-z ~ 2(1 -z) 



is factorized off explicitly and the factorization scale for 
is just Mi = a 2 of the gluon in the next K^f{k%). 

This fact that the factorization scale for the first emission is 
defined to be a 2 of the second emission is the essential differ- 
ence with the standard EGMPR/CFP scheme El|26l, where 
a, < [i for both emissions and therefore a 2 — > is not blocked 
by a i like here. The EGMPR arrangement has advantage of 
being similar to the system of UV subtractions |26 1 but causes 
oversub tractions, unfriendly for the MC implementation. We 
assume implicitly a cutoff regularizing a\ — > limit, for in- 
stance a i > «o. Note that, although the distribution of a\ 
seems to be simple, we cannot perform f =^ to get a pure 
log, because the upper limit a 2 > a\ is still nontrivial; and we 
have to wait until the next simplifications due to insertion of 
the second P', before getting the pure log from the integration. 

Let us now insert second P' into cf'V'K^V'K^). Again, 
we examine the limit a 2 — > 0, keeping 02/^1 = const. While 
taking this limit we keep a 2 = X\ — x 2 = const, such that 
$ = (P — k\ — k 2 ) 2 = const, in addition to previous s\ = const. 
More precisely, we start by introducing 



: s 2 



sx 2 = {P-k\ 



h-2(P-h)-k 2 



as an integration variable: 



f d 3 h d 3 k 2 , j 

o 2F = / dx\dx 2 —=^-— -Q-dx 2 [P -K\k\ -k 2 ;qi,q 2 ) 

Ci2>ct\ 



2C F a s P(l-ai) 

TZ 2 a\a 2 



pi (xipoF,PoB;k 2 ,qi,q 2 ) 



x 8^=1-01] s8(sx 2 -sx\ + 2(P-ki)-k 2 ). 



(22) 



Next, we perform the same transformations on 8-functions ac- 



' Here C is the interference term, while the other two terms are amplitude 
squares. 



companied by the rescaling k 2 = X 2 k 2 : 
s8(sx 2 — sxi + 2(P — ki) -k 2 ) 
= J dY s8(sx 2 — sx\ +2(P — ki) -k 2 ) Y~ 2 2poB-k 2 
x 8(5x2 — sx\ + 2pos ■ k 2 Y~ l ) 
dYs5(sx 2 -sxi + 2{P-k\)-Yk 2 ) Y~ l 2pwh. 
x 8(5x2 - ixi + 2poB ■ k 2 ) 



(23) 



^2 2poB ■ k 



x 2 -x\ 



2pm ■ k 2 



2{P-h)-k 2 \ s 
2Pob ■ h 



= Slx 2 -X\ + 

where 

X 2 (ki,k 2 ) = 



8*1 - *2=02J 

s(xi—x 2 ) _ s(xi_—x 2 ) 



2(P-ki)-k 2 2{P-X x {h)h)~k 2 



Note that the scaling factor X 2 — >• 1 in the collinear limit P2 — > 
0. 

Let us stress that the integral under consideration is now 
transformed into a new equivalent form, but the limit a 2 — > 
is yet to be taken! In the transformed variables the integral 
reads 



C2F = 



dxidx 2 d -^- d -^dz 2 {P-Xiki-X 2 k 2 -,q u q 2 ) 



ai>a\ >uq 

2C F a s P(l-a l ) 2 , . T , 

x — 5 ^— t — A-2 piixipoF ,P0B;Mk 2 ,qi,q 2 ) 

X ] ft, 8jq — X2=Oa ■ 



(24) 



Now we are ready to take the limit a 2 — > keeping a \ja 2 = 
const (also qq — > 0) and a, = const: 

C ¥ K 0F ¥ K QF - 

dxidx 2 -^--=^dx 2 (P-Xih-X 2 k 2 ',qi,q2) 

_ A j jLK^ 



M>a2>a\ 



2C F a s P(xi) 2C F a s P(x 2 xi) da B , 

X j ryy j -/? -377(^2,6) 

iz 2 a\a\ n 2 alai, dQ. 



2"2 



X 8_f 1 = l-ai 8^2=1— St— 82' 



(25) 



Note that the \\ factor from the phase space and the matrix 
element cancels out as before. The above formula is the prin- 
cipal result of this subsection. It defines the double use of 
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P', the transformation kj(kj) and its inverse kj(kj), i,j = 1,2. 
Note that in the above formulae we could use variables zi = 
1 — &i = x\ andz2 = (1 — &i — (Xi)/(1 — OCi) = X2/x\ instead 
of a,-, i = 1 , 2. In the following we may find it useful to switch 
to the Zi, i = 1,2, variables. 

With the global factorization scale M inserted at the 
end of the LO ladder the transverse-plane integration 



now decouples and provides a 



M>a2>a\>ao 

pure double lo^j 



C l?K QF lfK 0F 



In" 



M (2C F a s y 



2! «o V jt / 



J dx \Pqq O Pqq] 2R (x) <3 B (sx) , 



(26) 



where 



(c) A nice accident of the Jacobian \d(ki,k2)/d(ki,k2)\ being 
compensated by the matrix element is generally not guar- 
anteed. However, if it were not true, we would have to im- 
pose this by hand, such that a pure logarithm results from 



the phase space integration, like in eq. ( 26 1, similarly as in 
the definition of the collinear counterterm in refs. ll9l [3D . 
for example. 

(d) The role of the parametrization fc, — > kj in assuring uni- 
versality (process-independence) of the ladder parts will 
be clarified once the MC for the DY and DIS processes 
with the NLO corrections to the hard part are defined, see 
below. 

(e) The phase space parametrization in terms of kj will be 
used also inside Kq to parametrize the two-gluon phase 
space, for instace for the 2-gluon crossed diagram in the 
NLO ladder. 



4[<WJ° ) ] 2 s(z) 



1 



41n^+41n(l-z) 
o 



1-z 

-(l+z)lnz-2(l-z 



is just a double convolution of the LO kernel with the IR reg- 
ularization a,- > 8. 



The distribution of eq. ( 25 i is easy to generate in the Monte 
Carlo. First, one generates a, and a, paying attention to the 
constraint x — X2 — 1 — OCi — OL2, and are constructed in 
the laboratory frame. Then A,,- are calculated and the rescal- 
ing k^ — > k^ is done (in two steps!). Finally, in the frame 
P — k\—k2 one generates qt according to the Born differential 
distribution. The phase space boundary s > is obeyed auto- 
matically. In the MC the soft IR regulator z,- = Xijx^x < 1 — 8 
will be introduced and the overall virtual Sudakov form-factor 
?%fe ln M ln r 



exp( 



y will also be supplemented. 
The above example demonstrates the most important fea- 
tures of the P' projector (see ref. [16| for more details). In 
particular the following lessons are to be learnt: 

(a) The phase space parametrization ki — > ki plays an impor- 
tant role in P', as it is instrumental in implementing off- 
collinear extrapolation, and also helps to take the collinear 
limit in the first place. 

(b) The rescaled 4-momenta ki violate the 4-momentum con- 
servation, similarly like 4-momenta after action of the 
kinematical projector of refs. PI |27l . In our case, how- 
ever, the off-collinear extrapolation is effectively undoing 
this kinematical projection 11 ! an d allows to operate in the 
original phase space, with the 4-momentum conservation 
untouched. 



10 With our phase space parametrization in terms of ki the above mechanism 
of producing pure logs is a general phenomenon, because M is always at 
the end of the ladder, and the ladder has a built-in time-ordered exponential. 

" The undoing is "effective", because the kinematical projection of refs. |4; 
1271 in our methodology is never done. A clever parametrization of the 
phase space is the only thing really done. 



(f) The soft eikonal limit is protected by P', because rescal- 
ing ki ki preserves it. 

For a better (complete) understanding of the construction of 
P' one needs to examine in a fine detail the case of two gluons 
in the middle of the ladder, for instanc^j in CqP'Kqf (1 — 
W ! )Kqf, which provides the NLO correction to the evolution 
kernel. 



III. MONTE CARLO FOR EW BOSON PRODUCTION 

The insertion of the P' operator into the LO gluonstrahlung 
ladder with any number of gluons can be done similarly as 



in Sec. II B obtaining the distribution ready for the LO MC 
modelling of the production process of the EW boson with 
multiple gluons emitted from the incoming quarks. 



A. Simplified single ladder case 

We start with the gluonstrahlung ladder just in one hemi- 
sphere in order to avoid algebraic complications of the two- 
ladder case (to be dealt with in the next subsection): 



4 0) -r> = £ {cJVjsS?)"} 

77=1 



T.O. 



SF Hdx (fld^^J^Pizi) 

n=Q J \i=\ % 

n dOg 

x dx 2 (P~ *•/•;«!> #2) 9s<r| A =n » =1 , ; . ~^( sx >®)> 

(27) 



This task is pursued in separate works 1 1 311 1 VI - 
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where 



verse integration decouples and is feasible analytically: 



kj — hjkj, hi — 



*(*/-! -xi) 



2(P-r j r = \k j )-k i 

i i 
Xi = 1 - «i = Yl z J' Zi 



7=1 



Cq°' - Tf = J dx G F (M,x) Gb(sx), 



G F (M,x)=e 



It A « 



M 

9o 



(28) 



and P(z) = j(l +z 2 )- The a-ordering, M > a, > a,_i / 
I . ..../;. is rephrased into an equivalent rapidity ordering E < 
T| ( - < r|,- 1 < T|o, where M = y/se~~, see the definitions of the 
phase space-integration element c£"E(k) and of other kine- 
matic notations in Appendix [A] The T.O. subscript stands 
for the time-ordering exponential structure in the factorization 
scale, see ref. |[T6]l for a general definition. The S F function is 
the usual MC Sudakov form-factor depending on the shape of 
the IR-boundary a, — > 0. Factorization scale is now defined 
as the minimum rapidity E for gluons in the F hemisphere 
(maximum rapidity for gluons in the B hemisphere). For the 
moment S is a free parameter, to be defined more precisely 
later on. 

Note that the definition of A,,- is recursive, that is to define 
Xi one must know In a typical MC event first A,'s, cor- 

responding to very collinear gluons will be very close to one, 
X, ■ ~ 1 ; only the last ones corresponding to non-collinear non- 
soft gluons, i.e. close to the hard process, will be rescaled by 
a significant A,,- ^ 1 factor. 



where the IR regularization (1 — z, ) < A is used in the n-times 

(0)i®„ 



convolution of the LO kernel [P qq ] . It should be stressed 
that the above LO formula represents the LO MC without any 
approximation. 

The inclusive (bare) PDF G F (M,x) of the MC obeys by 
construction the LO evolution equation: 

fj G F (M,x) = [y$®G F (.M)](x), (29) 



3lnM 



1+z 2 



where ?f>(z) = 
the use of the T.O. exponent in eq. d27 



2(1-;) 



B. Two ladder LO case 



This is essentially due to 



Similarly as in the case of two gluons in eq. (26 1, the trans 



Let us now consider the case of two ladders. In the back- 
ward (B) hemisphere x, are related to p, (instead of a,), the 
evolution runs towards larger rapidity, otherwise all algebraic 
structure is the same. Again, we first express the LO MC mas- 
ter formula in terms of rescaled 4-momenta ki = A.,£,, but we 
postpone the definition of A,,, as it will be a little bit special. 
We propose the following multigluon distribution to be imple- 
mented in the LO approximationf^] 



n ] = 1 «2— 1 

«l=0 n 2 =0 J J^<r\, H \, =1 / 

f / "2 _ 1C F OL \ «l+«2 dOg 

Xe 5S X> (ll ' 32 ^) 9 ^^-!^^^)] h XB =li% { z B] d ^ P - E kj-q u q 2 ) ^{sx F x B ,Q). 

(30) 



In the above we understand that the F part is the forward part 
of the phase space T|of > T| ( - > E and the B part is the backward 
part E > T[i > T[ob- The rapidity boundary E between the F 



and B parts is kept as a free parameter as long as possible, 
to be adjusted later on. In particular E = is perfectly legal 
and may serve as the first choice, before something better is 
founcr^ Variables z F / and zbj could be defined similarly as in 



The reader should keep in mind that the above is for the "primordial" quark 
and antiquark initial beams, and their distributions in hadron will be added 
in the MC program. 



In practice the choice of S may be quite complicated, for instance it can 
be correlated with the rapidity position of the EW boson. In the NLO 
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the single hemisphere case, provided we perform kinematical 
mappings asymmetrically, first in one hemisphere and then 
in the other one. We shall do it below, however, in a more 
sophisticated symmetric way. The definition of zfi and zbj 
will result from that. 

The angle 8 can be defined with respect to any reasonable z- 
axis in the rest frame of P — for instance along pof — /%b 
in this frame. The boundary between the F and B phase space 
is at the rapidity S. It is also understood that the differential 
cross section of eq. < [30| ) is implicitly convoluted with some 
initial quark distributions at r\of and T|ob and the appropriate 
boost is done from the reference frame of the quark-antiquark 
frame pop + pos to the laboratory frame. 

Again, in eq. < |30) , phase space can be integrated analyti- 
cally over transverse momenta, providing the classical factor- 
ization formula 

C^T^r^ = / dx F dx B G F (Z,x F ) G b (Z,xb) C B (sx F x B ), 

(31) 

where G B (Z,x) — Gf(S,x). The remarkable feature is that 
the above LO formula represents the exact LO MC without 
any approximations ! 

We could define the dilatation parameters A./ (recursively) 
first for one ladder and then for the other, but this solution 
would be asymmetric. Instead, we define the dilatation pa- 
rameters in kt = Xki in a more sophisticated way, for both 
hemispheres simultaneously. For that purpose we introduce 
a new ordering (indexing) of gluons according to the distance 
|r|, — S| from the rapidity S of the EW bosorpj Formally, we 
define a permutation 

71 = {7Ii,7t2, ...,7C„ 1+ „ 2 } 

of all gluons in the F+B phase space, such that 

> -S|, i= l,...,m +n 2 , 

With the help of the above simultaneous "double ordering" in 
the F and B hemispheres, we define in a recursive way: 

km=hkTti, h = ^7~~pi , X 'l t j /= l,2,....«i +n 2 - 



2(P-Lyj l k n ,)-k 1li 



(32) 



where xi now is 



n 



(F.B)-Kj , 



where (F, B) means that we insert in the above product either 
ZFj or z B j, depending whether 7Cy points to the F or B region. 
The parameter A; is defined in eq. ( 32 1 recursively by means 
of solving the following equation 

si = sxi = (P-t k n f = (p-t M*y) 2 - 



The LO Monte Carlo algorithm using the above algebraic 
framework can be described step by step as follows: 

1. Variables x F and x B are generated with the help of the 
FOAM program [32 1, next parton (gluon) multiplicities 
rt\ 2 and variables ZFi and zbj are generated using the 
constrained MC (CMC) algorithm [33 10 The S vari- 
able is set. 

2. The four- momenta kf are defined separately in the F and 

B parts of the phase space using the CMC module, with 

the corresponding constraints: 1 — £ OLj = Y\ zfj—xf 

jeF jeF 

and 1 - I py = n z B] 



x B . 



3. Permutation n with the simultaneous ordering in F + B 
space is established. 

4. Using P and k ni the rescaling parameter Ai is calcu- 
lated, then^j =X\k lll is set. At this stage (P — k ni ) 2 = 
sxi , where x\ = z% 1 = 1 - oc^, or xi = z% x = 1 — Ptc, , de- 
pending whether k %l was in the F or B part of the phase 
space. 

5. Using P — k ni and k n2 the parameter A-2 is found and 
km = ^.2^2 is set > enforcing (P — k ni — k %1 ) 2 = sx2 — 
sz% l Zn 2 - The recursive procedure continues until the last 
gluon. 

6. In the rest frame of P = P — Y,jkx > lepton 4-momenta 
cfy and ^2 are generated according to the Born angular 
distribution. 

The definition of zf] and zbj in terms of Sty, j £ F and 
Py, j £ B, follows from the above algorithm and is more com- 
plicated that in the case of one hemisphere. The main advan- 
tage of the above scenario is that this way the kinematics of 
the two hemispheres is interrelated very gently, starting from 
very collinear gluons (for which A,,- ~ 1) and finishing with the 
least collinear ones, next to the hard process EW boson (the 
angular ordering is the same for kj and hi). 

As compared to ref. |[34|, where a similar algorithm based 
on CMC ll33l and rescaling of 4-momenta was proposed, the 
present algorithm does the "rescaling'f^] in a more sophisti- 
cated way. The rescaling affects mainly the hard non-collinear 
gluons, not the soft and collinear ones, while the rescaling in 
ref. Il34ll was global, similarly as global manipulations on the 
4-momenta (boosts and rescaling) used in other parton shower 
MCs flT] [8) in order to impose the 4-momentum conservation. 
Moreover, kinematic parametrization of the phase space in the 
present MC is based on the projector P' of the collinear factor- 
ization (instead of being ad hoc), which is essential for com- 
pleting the NLO. 



corrections calculation, for a single gluon, S = should be set. 
15 Which is also the phase space boundary between the F and B ladders. 



16 In the CMC algorithm parton multiplicity is generated as the first variable, 
contrary to the Markovian and backward evolution algorithms where parton 
multiplicity is generated at the end of the MC algorithm. 

17 The "rescaling" in our case is merely a synonym for parametrization of 
exact phase space in terms of the rescaled 4-momenta. 
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C. Real NLO correction to hard process 

The NLO correction in the EW boson production hard pro- 
cess (non-singlet) will be implemented using a single "mono- 



lithic" MC weight, which reweighs the LO distributions de- 
fined in the previous subsection. The real emission part of the 
NLO correction in this weight comes from the integrand of 
eq. ( 17 1 which we rewrite as follows: 



(l-p) 2 Ja B (l-a) 2 do B/ „ a 



'cop 



l + ( l-a-p) 2 da B 
2 d0. a 



(f,§)-e 



(33) 



l + (l-a-p) 2 da B 

~1> dlla 



(f,9). 



In the following use of function $\(p F ,p B ,...), defined in 
eq. (33 1, vectors pp and p B (p FB — 0) result from the last 
insertion of P' before the hard process and they are defined 
in the rest frame of P = q\ + qi to determine = Z(q\,po F ) 
in the LO part of the Born cross section. On the other hand, 
angles Q F and Q B are defined with respect to the original —poB 
and pop in this frame. They will all coincide when all gluons 
get collinear. 

The distribution of the MC with the LO+NLO hard process 
is now defined as follows: 



,(i )r (i) r (i) 



H / dxp dXB 

«1=0 n 2 =0 J 



jjd^m^,^,- 

<r h \i=i 



2C F a s 



"2 



2C F a s 



'*>r\n 2 V;=l 



where the MC weight is 



dc B {sx F x B ,Q) nrNW 



d£l 



W, 



MC > 



(34) 



NLO 



MC 



l+As+v+Y 



$\(s,PF,pB;aj,z F j) 



j^pP(zpj)dc B (s,Q)/d£l 

y &l($,pF,pB;aj,Z B j) 

£ B P{z Bj )d(5 B (§,Q)/d& 



(35) 



with the NLO soft+virtual correction Ay + $ to be defined sep- 
arately in the following subsection. 



Our construction of the above MC weight of eq. ( 35 1 is 



quite similar to that proposed in ref. [ 35 ] for the NLO correc- 
tions in the middle of the ladder^] The difference is that the 



18 The MC weight of ref. 1 35 ], when applied to the hard process, would also 
render complete NLO in the hard process and provide the same inclusive 
LO+NLO cross section. 



I 

proposal of ref. 1 35 1 was based entirely on the Bose-Einstein 
symmetrization of the multigluon emission in the LO MC and 
the resulting weight was more complicated, while the MC 



weight of eq. (35i is significantly simpler and algebraically 
similar to that of ref. [36| (albeit for the collinear rather than 
soft gluon resummation). 



The meaning of the arguments in pi is such that in eq. ( 33 1 
§ = sx F x B , and three angles 0, Q F and Q B have already been 
explained above. What remains to be specified is the defini- 
tion of a.j and p ; - in terms of the variables aj,Zj in j 6 F,B 
parts of the phase space as follows: 

CLj = l-zpj, Py = Cj aj /al , for j e F, 



$j = l-ZBj, 



0C; 



P 7 a|/a 2 , for j G B, 



where the rapidity r| ; is translated properly into aj and the 
rapidity E corresponds to as. 

As compared to earlier attempts to consistently implement 
the NLO corrections to the hard process in the parton shower 
MC, the proposals of refs. [34 , 37 1 were going in a similar di- 
rection. However, the present work differs from these works 
in three important points: (i) the virtual corrections A$+y are 
added here, (ii) the method of combining the NLO correc- 
tion with the LO MC proposed here is systematic and NLO- 
complete, (iii) the treatment of the kinematics is compatible 
with the principles of the collinear factorization. 

Comparison of our methodology with the well established 
met hods o f MC@NLO (9| and POWHEG flU is done in sec- 
tion nun 



D. NLO analytical factorization formula 



For the LO MC defined in eq. ( 30 1 we have seen that with- 
out compromising the exact phase space for multiple glu- 
ons (keeping the 4-momentum conservation) we could get the 
contributions of ladder parts factorizing off, exactly as in the 
traditional collinear factorization (in which the 4-momentum 
is not conserved). The above seems to be almost miraculous 
in view of the complicated nature of the exact complete phase 
space and the fact that no approximation was done. 
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What is even more amazing is that the same nice exact fac- 
torization is also true for the MC with the NLO-corrected hard 



process included according to eq. ( 34 1. The key point is that 



analytical integration of the phase space, again without any 



approximation, in eq. ( 34- 1 is feasible (!) and leads to a simple 
familiar result 



C«lf = J dx F dx B dz G F {Z.x F ) G B (Z,x B ) 



x o B (szx F x B )\ 8, = i(l + A s+V ) +C 2r (z) 



(36) 



where C 2r {z) = 2 ^ s - [-^(1 -z)] , see Appendix B. Two LO 
PDFs, G F (Z,x F ) and G B (Z,x B ) are these of eq. (31 1. The al- 
gebraic proof of the above formula can be found in Appendix 

<£E 

Note that the LO PDFs, G FB (Z,x), are by construction in 
the collinear factorization scheme specific for the MC, it will 
be the same in the MC for the DIS process - the scheme de- 
pendence of physical observables will cancel as usual, for in- 
stance while transferring experimental knowledge on the par- 
ton distributions from the DIS process to the DY process (or 
vice versa). 



Let us comment on the z-dependent C 2 , -(z) term in eq. ( 36 1, 
since it is different from what we see in ref. [300 where it 
is simply absent. It is not present there because it is com- 
pensated by the twin terms ~ — j£(l — z)l5(l — y 2 ) origi- 
nating from the y-traces and located at y = ±1, that is ex- 
actly at the collinear poles (for gluon strictly parallel with the 
quark). This cancellation is not disturbed by the MS collinear 
counterterm subtraction. In our MC factorization scheme this 
term is included in the counterterm, hence the net contribution 
from one of these terms, contrary to MS, stays uncancelled in 
eq. d36]l. 



under control: 

c£f M = [ l dzf$ M (z)o B (zs), 
Jo 

f^(z) = 8(l-z)+8(l-z)^ ( 2 -n 2 - 7 - 



C F a s ( S Y fP(z)\ fl 



K \fr- J \1— Z 



03 2 



(37) 



7C I 1 -Z 1 -Z 



where P(z) — \( \ +z 2 ), CO2 = Je — In 471 and the incoming 
quark and antiquark are massless and on-shell. The NLO 
real 1 -gluon emission distribution entering in the above AEM 
result is identical to Z}W of eq. (15i and 0^f M includes the 



complete all (gluonstrahlung) virtual corrections, including 
the quark wave-function renormalization constant Z F . Since 

°DY 



is gauge invariant, the calculation of eq. ( 37 1 is done in 



ref. Il30ll in the convenient Landau gaug^j 

In the formal standard MS methodology one subtracts from 
Idy M ( z ) two LO collinear counterterms 



(z) 



C F a s 



1 -zL e 



in order to avoid double counting with the ladder and/or ex- 
perimental PDF, thus obtaining Cq(1 — PA'of — VKq B ). This 
gives rise to the standard subtracted DY analog of the DIS 
coefficient function in the MS scheme 



Tdy (zj = Jdy (zj - K F (zj - K B (z 



= S(l-z) + 8(l-z) 



C F a, 



-7C 

3 4 



C F a s fP(z)\ /, § 

■ 2 ^r (i^J + (V +C02 



(39) 



■2P(z) 



lnz 

~z 



E. NLO soft+virtual corrections 



Let us consider the unsubtracted results of Altarelli-Ellis- 
Martinelli (AEM) work |3O0 in which real and virtual 1- 
gluon emission diagrams are combined and integrated over 
the phase space, keeping the variable z = x = s/s = 1 — a — p 



19 A quite similar formula for introducing the NLO correction to the LO 
kernel in the middle of the ladder was already tested numerically in the 
MC exercise with 4-digit precision, see ref. 1 16 1. 

20 It is also absent in QED, in the well known formula of Bonneau- 
Martin (38). 

21 See eq. (89) there. 



Note that the In -4 + co? term in the above will be absent if 

1 1 

the relation S = /j 2 e~ ai is adopted^] as in refs. ||39ll40ll . It is 
well known 113011 that the numerically dominant term '"/i^ in 
the above function is correcting for the misrepresenting of the 
soft gluon behaviour and incorrect phase space limits of the 
MS dimensional regularization (subtraction) scheme. In our 
MC scheme this term will not be present (gets transferred to 
the ladder). 

The plus regularization (...)+ of the IR singularity in 



eq. ( 38 1, in the diagrammatic approach of the CFP [26 1, comes 
from Z F in the axial gauge. It is also shown by CFP that there 
is a diagram-per-diagram correspondence between Feynman 



22 In Landau gauge Zp = 1 up to the NLO level. 

23 This corresponds to the use of 1 + CO2 instead of the pure pole in the coun- 
.,2 



terterm and setting s = fi 
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diagrams in the axial gauge and the diagrams of the operator 
product expansion (OPE) [41 , 42 F] 



In the context of the subtraction/factorization scheme 



aligned with the MC, we should apply eq. (17 1, that is in or- 



der to avoid double counting with the ladder (PDF), we should 



subtract from eq. ([6} the LO contribution of the MC: 

, 1 C F a s P(l-a-p) 
pi Fc (a,p) = s — i — % <a . 



p % 



a 



However, in order to combine it properly with the virtual 
corrections, the above distribution should be extrapolated to 
n = 4 + 2e and integrated over the phase space. This is done 
for the F hemisphere in the following (keeping in mind that 
1 - a- p = 1 - a): 



Cf a s 



CfO-s f dadfi 

"ajT 



2e 



sap 

Z/U 2 



P(l - a,e)0p <s 8a=i--e a>5 -5(1 -z) S M c{s,e) 



n 2+2E P(z,e) 



(1-z) 



1 r 2 i e i-z>8- 8 ( 1 -Z) Smc(s,£) = 



C F a s ( P'(z,e) 

K [ l-Z 



1 , S 

- + co 2 + In — r 
e z/j z 



(40) 



where the 0(e) contribution from y-trace is added to the LO 
kernel 



expanded Sudakov form-factor) of the LO MC reads: 



P(z,e) = P(z)(l+2eln(l -z)) + ^e(l -z) 2 , 
P(z,e)=P(z)+ 1 -e(l-z) 2 . 



The overall plus-prescription is coming from the 1st order ex- 
pansion of the Sudakov form-factor in the MC in n = 4 + 2e 
dimensions and from the usual sum rule J dz Kf c (z,e) = 
(treating S as z-independent): 



Smc (s,e) 



1-5 



dz 



n e Jo (1 -z) 1_2e \zfj 2 



(41) 

with both 0(e) terms in P'(z,e) necessarily participating in 
the (...)+ prescription^ 



The expression of eq. (37i is the result of the following 
phase space integration 

C F a s f dadfi f A«xp x £ 



fAEM/ P ^ C F a s fdadp r 

f D y (z,e) = — y ^-y^ 



2£ 1 9 

V 1 



Pi(a,P)5 1 _, =a+ ( 3 e a+ (5 >5 -8(l -z) U s+ v, 



(42) 



where pi (a, p) is the one real gluon distribution and Us+v 
sums up the real soft gluon, a + P > 8, and vertex virtual con- 
tribution. The complete 1 -gluon contribution (including the 



K^ c (z,e)+Kr(z,e) 



-mc/ Cf<x s f dadfi 



-dQ.2. 



2 c 



sap 



n J ap ' V ZjU* 
x [ J P(l-a,e)9p <a 8 1 _, = a +P(l - P,e)9 M 8 1 _ z= p ] 
x Si z>5 - 5(1 -z) 2S MC (s,e). 



(43) 



The difference between the complete NLO of eq. ( [37} and the 
above LO MC contribution, after partial phase space integra- 
tion, reads: 



Tdy (z) = Joy (z, e) - K F (z, e) - K B (z, e) 



CfV-s 



(l-z)+S(l-z)A v+ s, 



C F a, 2 2 7\ C F a s 1 C F a s f 2 , 5 

Ay +5 = 1 -K H 1 - = 1 -71 

+i 7C 1 3 4 i 7C 2 Jt 1 3 4 



(44) 



From the above we are able to determine the z-independent 
soft+virtual correction Av+s in the NLO MC weighPj The 
above does not include any singular terms like ln(l — z)/(l — 
z)+, as advertised earlier. 

The difference between the standard MS function of 



eq. p9\ and that of eq. (44 1 is entirely due to the difference 



between the MS counterterm of eq. ( 38 1 and the MC countert 



In this sense the axial gauge is implicitly present in eq. \3Sf , while the 
unsubtracted eq. |37| is gauge invariant. 

In the MC practice regularization of 1/(1— z) is done with some cut-off 
1 — z > § rather than with £ of the dimensional regularization. 



The last term in Ay+s is due to the plus-prescription in the (1 — z)+ part of 
the MC counterterm of eq. J40|. 
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erm of eq. (40i: 



f& S (z)-f&(z) = -2KP*(z,e) + 2KP(z,e) 



cMC 



MSi 



rMCt 



C F a,s , , ~C F a s ( P{z) 
(1-ZJ++2— — 

71 V 1 —Z 



7C 

C F a 



In -r + co 2 



r 4^) ln ^ ) -2P (z) i^ 

71 [ 1 — Z 1 — Z 



(45) 



and it represents clearly the difference between the MS and 
the MC factorization schemes. 

One may ask how to interpret this change from MS fac- 
torization scheme to MC factorization scheme, in particular, 
how unique the modified MC counterterms of eq. pO} are. 
One may answer this question in two complementary ways. 



One way is that the new MC counterterm of eq. (40 » repre- 
sents just the collinear limit of the exact matrix element in 
n = 4 + 2e dimensions (keeping higher order terms in e) in a 
sense of the P' projection operator. This definition has to be 
supplemented with the plus prescription in the soft limit or, 
alternatively, by saying that Zp, which in CFP (MS) provides 
for plus prescription, is replaced by the Sudakov form-factor. 
This approach represents an effort in combining the best from 
the two, the collinear and soft resummation. Another way of 
addressing this question is to say that the real backbone in 
the collinear factorization is OPE with CFP providing a solid 
bridge to OPE, and the only thing that has to be explained 
and kept track of is the difference between CFP and MC (in 
a similar way like finite UV renormalization). This approach 
was already advocated in refs. [4, 43 1 and in other papers [39], 
where factorization scheme dependence was discussed. In our 
approach, we are using both ways of addressing the above 
question. 



The related question is whether the counterterm of eq. (40 1 



is universal? Basically the answer is that it is universal thanks 
to the fact that it is defined in terms of the kf = Xk? variables. 
In other words, the kinematic mapping, inherent in the new 
P' operator, should remove the hard process dependence on 
the side of the ladder, in the same way as the pole-part oper- 
ation in CFP ll26l or P^m of ref. (4) . To be completely certain 
that the above aim of the universality of the new MC factor- 
ization scheme is achieved, in next section we shall define a 
similar MC scheme for the DIS process, define and use the 
collinear counterterm of this MC scheme, and in Sect. |IV I| 
we shall check validity of the factorization scheme indepen- 
dent relation (DY-2xDIS) of ref. (30) between the coeffi- 
cient functions of DY and DIS, both taken in the MC factor- 
ization scheme. 



F. Differences compared with POWHEG and MC @NLO 
methods 

In this subsection we outline main differences of our 
method compare to the well established approaches of 
POWHEG [11| and MC@NLO used today to combine 
NLO-corrected hard process with the LO parton shower. 



First and most obvious difference between our method 
and those of POWHEG and MC@NLO is use of different 
factorization schemes. In our approach we use factoriza- 
tion scheme |[T6l [TTIl designed especially for MC simula- 
tions whereas POWHEG and MC@NLO use the standard 
MS scheme. This allows them to use directly the standard 
MS collinear PDFs, while we need additional work herq^l 
Moreover, we build the LO parton shower MC from scratch 
whereas POWHEG and MC@NLO profits from the well es- 
tablished (unmodified) LO MC programs. 

At first it may seem that these general features result in un- 
necessary complications in our approach, however, profits are 
more important, especially if we have in mind construction of 
the fully NLO MC parton shower (with NLO corrections, not 
only in the hard process, but also in the ladder parts). Our 
method features: 

(i) simple and positive MC weight implementing the NLO 
on top of the LO MC, see ref. ED (MC@NLO features 
negative weights); 

(ii) no need to correct for the difference in the collinear 
counterterm of the LO MC and the standard MS scheme; 

(iii) the virtual+soft corrections Ay+s are completely kine- 
matics independent - all annoying dL c± contributions of 
MC@NLO are gone; 

In" ( 1 — x) 

(iv) built-in resummation of the — j _ ' terms is provided; 

(v) direct relation to the collinear factorization procedures. 

Note also that, in the presented method there is no need to 
define the hardest emission, as in POWHEG, as it is automat- 
ically included into the sum over spectator gluons in the for- 
mula for the MC weight in eq. ( |35] >. In fact, we can explicitly 
see that the dominant contribution is from the "hardest" (in 
k T ) gluoiQ for numerical illustration see ref. [21 1. This al- 
lows us to avoid truncated/vetoed gluons needed in POWHEG 
methodology in case of angular ordering. 

A detailed comparison of MC@NLO and POWHEG meth- 
ods themselves can be found in refs. fi4l l45ll . 



IV. DEEP INELASTIC ELECTRON-PROTON 
SCATTERING 

As already said, the process of deep inelastic electron- 
proton scattering (DIS) is included in the scope because it is 
is an important source of information on parton distributions 
in a proton and, by comparing the DIS and DY processes, the 
question of universality in the collinear factorization imple- 
mented in the MC can be fully discussed. 



27 One possibility is refitting PDFs, which shouldn't be to complicated as the 
difference between MC and MS scheme on the inclusive level is small. 
This is just a relabeling according to kj, we do not need to change previ- 
ously generated, angular-ordered gluons. It is only exploited here for the 
purpose of efficient evaluation of the NLO MC weight. 



2S 



16 



In the following subsections we shall first introduce kine- 
matics, phase space and notation for one real gluon emission. 
Next, we shall define the multigluon LO MC distribution with 
initial-state radiation (ISR) and final-state radiation (FSR) LO 
ladders and the LO matrix element for the hard process for 
electron-hadron DIS. Analytical integration of the MC distri- 
bution will lead to the familiar formula for structure function 
F2 in the form of the convolution of PDF with the Born cross 
section. Then we shall give a close simple formula for the 
MC weight implementing the NLO correction to the hard pro- 
cess. The analytical integration will be again possible giving 
the structure functions F2 and F\ in form of the convolution 
of PDF with the NLO coefficient function. Of course, the 
above NLO coefficient function will be in the MC factoriza- 
tion scheme, but we shall see that universality is preserved, by 
means of checking validity of the factorization scheme inde- 
pendent relation DY— 2xDIS of ref. [ 30 1 between the coeffi- 
cient functions of the DY and DIS processes. 



A. One-real gluon distribution and kinematics 







U=q 2 +kj 




Pi 








• cm y 












P 2 













The differential distribution in this case reads 



Qq a QED 



s z + u 



,s/ 2ffi 
C F a s dadfi t\ 
71 ap t ' 



d\\f 
2% 



(47) 



The additional invariants s\ = 2p2qi, «i = 2q\p2, ?i = 2*71^2 
are introduced in this case. The factor ~ is the Jacobian due 
to parametrization of the phase space in terms of the rescaled 
Sudakov variables 114611471 : 



a : 



2kq2 
t\ +2kq\ 



2kqi 
t\ + 2kqi 



(48) 



The angle \|/ is the azimuthal angle of k around z-axis in the 
Breit frame of Q — q2 + k — q\ , that is where Q° — 0, with an 
additional requirement that q\ is parallel to the z axis. We call 
this the reference frame B, see Fig.|4] 

Yet another Breit frame B\ is marked in Fig.[4j that of Q\ = 
qi — qi with the z-axis along q\ . It will be used in the MC and 
in the analytical calculations. Note that the integration is over 
the angle \\t of the k,qi,q2 plane as a whole around fl, while 
another azimuthal angle (j)i of k in the B\ frame is frozen to 
zero. Note that the standard Sudakov variables are: 



kq 2 



kq i 

qiqi qiqi 



a' = ^, p' = ^-, t = fi(l-a' + p')- 



(49) 



They are not convenient because p' e (0, °°), and the following 



transformation is mandator 
a' „ p' 



a = 



1 + P'' 



1 + P' : 



a' = 



a 



P' = 



P 



0<oc<l--, < P < 1, 
s 



1-P' ' 1-p' 



(50) 



B. Bjorken variables, structure functions, collinear limits 



FIG. 4. Kinematics of one gluon emission in the Breit frame. 

The Born differential cross section of the electron-quark 
scatterning e(p_A + q(q\) — > e(p2) + q(q2) in terms of the stan- 
dard variables; | s = 2p\q\, t = 2p\p2, u — 2p\q2 reads 



. r , 
dCs = — d 

s \s 



s 2 + u 2 



(46) 



where Q q is the quark charge. Next, consider the process with 
emission of an additional gluon from the quark line: 

e{pi)+q{q\) -^e(p 2 )+q(q2)+g(k). 



The standard Bjorken variables are 

t ie 21 



2q\Q 2q\Q 



, 1 > x B > 0, y B 



Qqi 

p\q\ sx B 



(51) 



In the case of a single gluon they are expressed as follows 

t l-oc' + p' , 2^ig t 

x b = ^ — ^ = — r - 7T, — = 1 - a, y B ~ - 



2q Y Q 1 + P' 



2q\p\ s(l-a) 
(52) 

The reader should keep in mind that for simplicity jcg is the 
fraction of the parton momentum in the initial quarl 



yxn 
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We omit the minus sign in the variables like t and u with respect to the 
standard notation. 



In the collinear limit k~ q2$/(\ — $), with z = 1 — P being the lightcone 
variable in the LO splitting kernel. 

Returning to the normal definition in the MC (a fraction of the hadron mo- 
mentum) is quite trivial. 
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Let us recall the definitions of the standard deep-inelastic 
structure functions in terms of the above Bjorken variables: 



C. Bare structure functions for 1-gluon emission 



d 2 a _2%v 2 QED Q 2 q 



dtdxg 



{y 2 B 2x B Fi(x B ) +2(1 -y B )F2(x B )} , 



2na 2 om Q, 



QED^-q _i 



27ia; 



*B l {I 1 + (! -ysf] F2(x B )-y 2 B x B F L (x B )} , 



+ 2(1 -jbW^l^b)}, 



(53) 



where we have employed the standard definition 2xF\ = Fi — 
xFl. In the LO case the Callan-Gross relation 2xF\ = Fj is 
fulfilled and the longitudinal structure function Fi = (it will 
receive a nonzero contribution at NLO). The LO relation to the 
parton distribution function (luminosity) is 2F\(x) = F 2 /x = 
PDF(x). 

It is instructive to investigate the collinear ISR and FSR 
limits. The slightly reorganized single-gluon emission distri- 
bution reads 



Our immediate aim is now to reproduce the well known 11261 
l48l l49l result for the NLO correction to the F2 (x) structure 
function by means of integration of the one-gluon phase space 
(the NLO correction to Fl(x) will also be found). The aim is 
to test our Monte Carlo phase space parametrization, prepare 
ground for determining the soft+virtual correction in the MC 
and put FSR under control. 

The unsubtracted (bare) contribution to F2(x)/x corrected 
due to the real gluon emissions plus the vertex correction 
(eq. (59) in 1 30 ] ) can be rewritten as: 



C F a s 

K 

f(l-z) 
z/n 2 



P qq (z) 



rl 

£ 



C0 2 



(58) 



, , dt d\\t C F oc s dadB 
d^ = a 2 QED Q 2 lT2 d^— -^w, 



W = 



s 2 + u\ + s\ + u 2 



(54) 



2s 2 



The soft limit is already manifest in the eikonal phase space 
factor . The following explicit expressions for the invari- 
ants in terms of our Sudakov variables are useful: 



where P qq {z) = 2(i~- z ) ' 032 ~ ^ E ~ m (4ft) comes from the 
&2+2z = 2tc(1 + £0)2 + • • • ) expansion. The baryon number 
conservation sum rule Jq dz C^We (z) = 1 holds explicitly. 

The standard NLO MS correction C\ to z~ 1 F2(z) form- 
factor is obtained simply by means of subtracting the MS 
collinear counterterm 32 i{P w (z)} + (i.e. the pole part). The 
formula of ref. |48 1 tobe reproduced reads: 



- = (i-P)y«, - = i-y B , 

s s 
5.~(l-a)(l-P)+aP(l-y B ) 



^-(l-a-PJy* 
s s 



AF 2 NLO (x B ) = C 2 (x B ) 



2cos ¥v /(l-«)(l-P) ap(l-y fi ) 



(55) 



In the FSR collinear limit, a ~ 0, p ~ 1 - z, k ~ qi ( 1 - z) /z 
and y B — yo — t/s, we have 

s 2 + u\ ~s 2 + {s-t) 2 , s\ + u 2 ~ (s 2 + {s - 1 f) (1 - p) 2 , 



W : 



l + (l-p) 2 ^ 2 + (^-f 



[l + (l-yo) 2 



(56) 



In the ISR collinear limit P ~ 0, a ~ 1 — z and k ~ (1 — z)q\ 
we have: 



1-a 



s 2 + u 2 ~[s{\-a)] 2 + [s(\-a)-t] 2 
l + (l-a) 2 r , n9i 1 



[i + (i-y») 2 ]- 

(57) 



p qq ( x B ) 
p qq {x B ) 



C F a s 

C F a s 

C F oC v f 1 — x B 
jt 



In 



In 



r(l-^fl) 

/J 2 X B 

t(l-x B ) 

fJ 2 X B 



0)2 



co 2 



3 1 

4 1-V 

3 1 

4 1-x B 



x B 
3 



(59) 



The last term (1 — x B )/2 in the non-singular part is due to the 
£-term from the y-trace for the intital-state collinear singular- 
ity and CO2 from the n-dimensional phase space. 

We start from the unsubtracted (bare) DIS distribution com- 
ing from two real-gluon emission diagrams from the quark 
line plus the vertex virtual correction in n = 4 + 2e dimen- 



32 We subtract a pure pole as in the original CFP work and not ( 5 + 
Gte){Pqq{z)}+, as it is a common practice nowadays. 
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In (1-x) 



soft 




► Ina 



FIG. 5. The logarithmic Sudakov plane for 1-gluon emission in the 
DIS process. 



sions: 

a 



d 2 a 2nan Fn Qf, r 



G 2 {xB,yii) = J dadfi J "" 2 2 + 2e p 2 (a,p) 8 tB= i_a, 
P2(a ' p) = ^T^ W(a ' P '^' £) U-a)(l-P) 

x 9l>a>59l>p 

s 2 + u\ + s\ + u 2 s 2 + u\(t-h) 2 



W(a,p,y B ,e) = - ^- 



(60) 



where we have reinstalled in W the £ term from the y-trac^] 
The real emission phase space is explicitly integrated for 
a > 8, where 8 -C 1 is an IR cut-off. The above phase space 
division is graphically shown in Fig. [5] Note that the a > 8 
part of the phase space, which we are going to integrate over 
includes not only the hard collinear ISR but also the hard 
collinear FSR! The costant U must include the vertex cor- 
rection summed up with the soft real emission a < 8. For 
determining U it will be enough to know |48 49 1 that the F 2 



part of the distribution in eq. (6O1 fulfills exactly the Adler 



we may omit the details of its calculation. The complicated 
phase space factor is simply due to the fact that the transverse 
momentum of the gluon in the Breit frame is: 



rap 



(l-a)(l-P)' 

Finally, / dQ^,^ is the n-dimensional extension of J ( 2jl d\\f. 

In the CFP scheme the ISR collinear singularity upon inte- 
gration gives rise to the LO pole part 



tfpAb, = cW fcO = |[i+(i -ys) 2 } ^ ( P(xs) 



1 -xb 



where Kqi is the lowest-order 2PI kernel for the ISR ladder 
and the + prescription comes from Z/r, as usual. The sub- 
tracted hard process matrix element in the CFP scheme is 
C o ] ~ C§ ] ~PKo,. We shall calculate it with the help of the 
usual counterterm technique. The ISR collinear/soft countert- 
erm (SCC) we define as follows: 

2l C F a s P(l-a) 



p 2t (a,p) = [l + (l-y B ) 2 ] 

X p £_1 B _e 9p<B(a)9l>a>5 5 

^(1-a) 



(61) 



B(a) = 



fa 



It is defined such that it integrates to the pure pole part exactly: 
cftf(x B ) = J </a£/pp 2c (a,P)6i- fl , =0 . 

In Fig. |5] we have also marked the integration area for the 
above counterterm. As we see, in this area the upper phase 
space integration limit from the energy momentum conserva- 
tion is replaced by the limit on the gluon transverse momen- 
tum equal (approximately) to fj. 

With the help of the above ISR collinear counterterm our 
task is reduced to calculate the subtracted DIS distribution in 
n = 4 dimensions: 



2n 



d\\i 



sum rule in the dimensional regularization, and in this way 



G 2 (x B ,yB) = J dadfi J — [p 2 (a,p)-p 2 c(a ! p)]8 TS= i_ a , 

(62) 

except for the trivial 8 term in W, which contributes 

^ir 1 (^T a )+- to be added at the end - The same with the 
similar ~ e<X)2 term from the phase space. 

The integration can be summarized as follows: 



1 Only for the ISR collinear singularity; the one for FSR falls into the U 



factor. 
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d 2 ^,° 2™ 2 qed% 2C F a s , d3 



dtdxB 

CI(xb) 
C l {xb) 



W(a, hy B ) - W (y B )P(l - a)0 Pl<B 



(«i) 



2na 2 QED Q 2 



([l + (l-y B ) 2 ]C s 2 (x B )-y 2 C L (x B )), 



C F a s 

K 

C F a s 

K 



P qq{ x B) 



In 



f(l -x B ) 

jA 2 X B 



CO? 



-x B , 



3 1 ,3 

Zi + l + o XB 

4 1-x B 2 



W Q (y) = l + (l-y) 2 , 



8,_ 



l-x B =a\ 



C F a s f 1 — x B 



(63) 



where the + prescription is provided by the virtual correc- 
tions. We have also included the e-contribution from the y- 
trace and ©2 from the phase space. As we see, C 2 (z) is equal 



from the different cut-off in the ISR counterterms: 



to the finite part of C^u^Jz) of eq. (58 I, thus we have repro- 



duced the classic result 11481 . as promised. 

In the MC scheme the ISR counterterm Cq^P'-Ko/ is defined 
as the single-gluon distribution which is extrapolated to n = 
4 + 2e dimensions and integrated over the phase space: 



C F a s rdadfi r ( rap 

Ki(z,£) = / — ^ / d£l 2+2 A T, r-» 

it i aP i V(l-a)^ 2 

x P(l - a,e) p<a 5i_ z=a 9 a>5 - 5 e= i5/ 
n 2 +2 g P'(z,e) 

u^9i-z>8 - 5(1 - z) S, 



(64) 



where P'(z,e) = P(z) (l +2eln(l -z)) + |e(l -z) 2 . The 
source of the + prescription in this case is the MC Sudakov 
form-factor calculated in n = 4 + 2e in such a way that the 
sum rule / dzKj(z,E) — is preserved also in n-dimensions: 




j*(z,e) 



'~ n £ Jo ^(l-z) 1 -^ 



hence two 0(e) terms in P'(z,e) necessarily participate in the 
(...)+ prescription. 



Subtracting Ki(z,e) of eq. (64 1 from the complete O(a') 
result of eq. (58 1 gives us the following coefficient function 



Cf C (z) 



C F a s 

K 



1+z 2 . ( . , 3 

■— rlnll — z) 

2(1 -z) V 7 41 



(65) 

in the MC factorization scheme, with the angular ordering. 

The most important part of the difference between the 
above MC structure functions and the MS variant is coming 



2%a 2 }ED Q 2 , c F a s 



d^ik^WoM 



-P(l-a)e Pl<B(ai) +P(l-a)e Pl<ai 



n-x B =o.\ 



2 ^ 2 qedQ] 



W (y B ) 



C F a s 



2P qq {x B )ln(l -x B )+P qq {x B )\a 



fJ 2 X B 



(66) 



A few comments are in order: 



Why not k -ordering? In such a DIS-like factorization 

would have been gone from 



scheme 



eq. (65 1 
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term 



ln(l-x 

l-x 



FSR could be treated in the DIS MC with- 
out the LO resummation, with the unexponentiated FSR 
NLO corrections. However, if the universality is to be 
maintained, and the same k~ -ordering is applied to the 
W/Z production process, that would either mean asym- 
metric treatment of the emission from the quark and an- 
tiquark lines or a large double logarithmic dead zone in 
the corresponding LO MC, between the ISR and FSR 
phase spaces. Both options are unacceptable. 

• Is there also a kinematic mapping involved in the above 
P', like in the previous W/Z production process? Yes, 
it is implicitly included in the definition of a and p vari- 
ables in eq. ( 50 1, where dilatation using factor 1/(1 + 
P') factor is seen. 

• From the point of view of the MC the above consid- 
erations are incomplete, as they still keep FSR in the 
inclusive/integrated form. 



D. DIS multigluon LO Monte Carlo 

Let us start with the raw distribution for n gluons, the 
(aC F ) n part only, relevant for the LO MC. 

e(Pl)+q{q\) -> e(p 2 )+q{q2) +g(k\)+g(k 2 ) + ■■■ +g(k n ). 



34 In the DIS factorization scheme C2 = exactly, while in the k ordering it 
would be only less singular. 
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The corresponding differential distribution reads 

W dy fC F a s \" 



do n = Q 2 q a 2 QED dtd<$ 



th 2% \ % 



(67) 



where W is a mild function to be defined later on. The in- 
variants s\ = 2/?2<?2> «i = 2</i/?2> h = 1q\q2 are the same as 
previously. The factor j is again the Jacobian due to the 
parametrization of the phase space in terms of the Sudakov 
variables ll50l . see bellow. The angle \\t is the azimuthal angle 
of k around the z-axis in the Breit frame of Q = qi + k — q\ 
with Q° = 0, with the additional requirement that q\ is paral- 
lel to the z axis. We call this reference frame B. Another Breit 
frame B\ is used in the MC, that of Q\ = q 2 — q\ with z-axis 
also along q\ (and qi). The illustration of the kinematics in 
Fig.pflis still valid, provided we replace k\ by 

The integration is done over the angle \|/ of the (II, #1,^2) 
plane as a whole around II = qi + while there is a sin- 
gle restriction on n azimuthal angles (]), of kj in the B\ frame, 
namely the vector YJcj must be co-planar with p\ and p 2 . 

The standard Sudakov variables are: 



, _ kicn , _ hqi 

u-i ~ ; Pi ; 

qiq2 qiqi 

j j j 

Next, we transform them as follows 1 50 1 



(68) 



a 



a 



t 
s 



(69) 



0<£a,<l- ? 0<£P ; <1. 



The Bjorken variable xg (of parton in the initial quark) can 
be expressed in terms of the Sudakov variables. Using Q = 
K + q2~qi and K = J\- kj, we obtain 

= 2q l q2+2qyK-2q 2 K-K 2 = ' ' I ; |3', I;^ „ 2 
2q 1 q 2 +2q 1 K ~ 



= 1 -Va r f = — — - — -j s 



(70) 



In the NLO world the term 

- 2 K 2 
K = 



2qiQ 



2 Y*i> j H ' kj 
2qxQ 



can be either omitted or taken care of in the collinear limit. 
Note that K 2 is absent in the case of the single-gluon calcula- 
tion of the NLO coefficient function. 



The fully differential distribution for emitting n gluons in 
the LO MC for DIS we define as follows: 



1 _ s dy 

Q ED uimy — 



da„ = Q 2 <X 2 QED dtdy \ e 



/ " C F a s dUidfr dtyi \ 



(71) 



The key objects to be defined are the variables 2j and the 
Sudakov form-factor S. For this LO modeling of the gluon- 
strahlung we use ordering according to the factorization scale 
(evolution) variable 



a ; V A 



(72) 



which is the variable of the angular ordering of the MC. 

The ISR part of the Sudakov plane (the blue trapezoid 
in Fig. |6| contains the gluons / = (1,2,3, ...m) which have 
P,/a,- < e~ and the FSR part (the red trapezoid in Fig.[6]l hosts 
the gluons f = (m + l,m + 2, ...,«) which have pi/a,- > e a . 
We shall indicate that the gluon j belongs to one of these two 
subsets by j E I or j 6 f . The variable z,j of the ISR and FSR 
gluon is defined in terms of either a's or p's: 



for jel: Zj^z'j = xj= l-£ 



a, 



;=1 



(73) 



for jef: ZJ = z] = 



K j-1 



1- E P« 

(=»»+ 1 



The Sudakov form-factor S is the integral over the area in the 
logarithmic Sudakov plane available for the real emission in 
the step-by-step Markovian process. This area is visualized 
in Fig. [6] as a shaded polygon. The rapidity S defines the 
boundary between the ISR and FSR emissions according to 
the corresponding LO distribution and can be treated as an 
arbitrary parameter, for example S = in Fig. [6^ is an ac- 
ceptable LO choice. In fact, another more clever choice of S, 
like the one indicated in Fig. ^jp, can be made, for instance, 
within the Markovian LO MC algorithm. It is also possible to 
switch from one value of E to another in the final stage of the 
LO MC by means of reweighting MC events. In view of the 
above flexibility, we leave out the exact definition of E to the 
later stage of the MC code implementation. 



E. Structure function for LO MC 

The standard double-differential distribution of the DIS 
process, as realized in our LO MC, is obtained by inserting 



21 




FIG. 6. The Sudakov plane of the LO MC for DIS. The shaded area denotes the integration domain for the Sudakov form-factor 5 of the LO 
MC. Rapidity equal to S marks the boundary between ISR and FSR. 



the 8-function defining xb 

d^ Lo _ 2*Qy QED 



dtdxs t 



E fl^f^^w. (74) 



n=0 V=l 

X 6 S 8 E^ : =0 e E«i<l e EP,<l 5 ^=l-Ea J -^2- 
The above distribution is directly implementable in the MC 



form, for instance using the Markovian algorithm, and the 
above double differential distribution of xb and t is com- 
ing just from histogramming, using MC events with all 4- 
momenta of all leptons, quarks and gluons explicitly defined. 

On the other hand, we may explicitly show analytically that 
the above distribution is proportional to PDF convoluted with 
the coefficient function. This is done by means of inserting 
into the integrand 1 =n?=i {®aj>a- + 0a,-<«*) an d after expand- 
ing/reordering the sums of the integrals. The distribution of 



eq. (74 1 almost factorizes into the ISR and FSR parts: 



d 2 a _2%Q-a 2 QED 



dtdxs 



in? 



„=o v=i 



in 

n'=0 \i=l 



2Cfa s f 3 



d "E(ki) P(zf ) a ,- 



d 3 <E{ki)P{zl)Qa s > ai >a^ 



EP iE 5<i-P/ 



Si 



9 E aj<A B =l- E aj 



(75) 



where we have split the Sudakov form-factor into the ISR and 
FSR parts, S = Si + Sf and P/ = £P, G /. In the above we 
have also neglected K 2 in xb — 1 — Y,j a j ~ K 2 . This is well 
justified at LO, but it turns out that it can be done at NLO as 
well. The alternative solution would be to make special effort 
in parametrizing the phase space (part of the definition of the 
P' operator) to "protect" xb as it was done for x = s/s in the 



W/Z production process. We have decided that this is not 
worth the effort, as the dependence on xb of the differential 
distributions is relatively mild. We may come back to this 
idea if an additional justification is found. 

Altogether, the final LO formula can be written as a convo- 
lution of the PDF for ISR and the the resummed "coefficient 
function" C<f(zf) for FSR: 
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d 2 a 2%Q 2 / a 



QED 



dtdXB 



Wo(y B ) J dxidz F 8 XB=XlZp D I (E,xi)Cjr(z F ), 



r I n 2Cf(Xs 



D I (E,x I )=e- s 'l i n 

n=0 \/=i 

{ „'=1 \i=l n J 

[ dz F C 7 {z F ) 
Jo 



d 3 <E(k) P (z/)e fl3 > a ,. >fli _ ] 8*=i_ Ey6 , 0y , 



I )^r)e« i > ai - 1 > as V 



(76) 



The interesting pure FSR object Cf{x) is probing the FSR 
evolution variable, instead of the FSR lightcone variable! In 
the LO version it is enough to keep only the trivial Cj (x) — 
8(1 — x) term, while for our NLO purpose it is enough to retain 
only one more easily calculable term, n' = 1, S = 0: 



d 1) (x) = 8(i-x) 



C F a 



Jo -sipr^-po^si-*^ 



8(1 



C F a s 



- ln( ^ ) - 3 ^ + 1 (3 + x; 

1 — jc 4 1-jc 4 V ' 



(77) 



NB. the above reproduces the bulk of the coefficient function 



of eq. (65 1, that is terms like ( '"j 1 ^ ) and 



The MC initial-state PDF obeys the LO DGLAP evolution 
equation (limited to the non-singlet gluonstrahlung): 

2tj t D { (t,x) = ^D { (t,x) 



I 



C F a s (P{z) \ 
dzdxb Xl=xz Di{t,x), 



l-z 



(78) 



and the same is true for the structure function 2F\ = C<y ®Di\ 



2t^F l (t,x) = ^F l (t,x) 
— / dzdx 8v, = 



C F a s f P(z) 
l-z 



(79) 



Fi(t,x). 



F. Exclusive ISR and FSR subtractions in DIS 

The two soft counterterms, for ISR and FSR, can be iden- 
tified in the fully differential distribution of the single-real 
gluon in the LO MC: 



, 1 dMf / 7 \ C F a s daid$i dtyi 

W (y B ) {P(l - ai)8 Pl<ai +P(l - pi)8 Pl>ai } , 



2% 



(80) 



where we define the y# = ,* * variable and the Born spin 

factor is Wo(y) = 1 + (1 -y) 2 . On the other hand, the NLO- 
complete unsubtracted distribution is 



do N x LO = Q 2 q a 2 QED dtdis? I g 8,(1 



C F a s daid$\ d§\ 



W(ai,Pi )yi ,) = 



7i aipi 2% 

s 2 +u\+s\ + u 2 
2s 2 



(81) 



See eq. ( f55] > for explicit Mandelstam invariants. 

For the MC we shall use the subtracted distribution with 
both the ISR and FSR counterterms: 



daf NLO = da^ w -daf cw 
1 dy 



C F a s da x d$ x <% 5 , 
i j — — PiW, 



= Q- q a QED dtd^ Tl -%\K j n ^ 
PlW = P/6p 1 <a 1 +Mp 1 >a 1 , 

P/(ai , Pi,y a ) = W(ai, Pi,y B ) - W (>> B )P(1 - oti), 
p F (a 1 ,p 1 ,y B ) = W(a li p 1 ,y B )-Wo(y B )P(l-Pi), 

which defines (up to NLO) the following expression 
~ C F a 4 A 1 A 1 rfai^Pi 

^ANLO(Z,yB/Z) — 



(82) 



71 Jo Jo aipi 27c 

x {P/Qp^a, +Pf9p 1 >a 1 } 8i- z=a , 

;r(l+z)ln(l -z) + -z+- 
2 4 4 



(83) 



to be used in the numerical tests of the MC implementations. 



G. Exclusive NLO correction to hard process in DIS MC 

In the following we propose a MC weight which upgrades 
the MC with the LO hard process and the LO evolution kernels 
to the MC with the NLO hard process and the LO evolution 
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kernels. The distribution in the LO+NLO MC reads 



where 



do N „ 



LO 



dt 



d\\i 
2% 



"C F a s daid$,d<bi- k \ 

xe E « i <ie E p. <ie - s w ()'B) w^ c LO , 

where the key element is the following MC weight 



(84) 



< c i0 = [l+A s+y ] + £ 



/e 5 W (y B )P( Zj -) + ; |^ Wo(yj)^(z;) 
-ti.ua i i V ^ a h z j^ i y M a hzJ >yz) 

- [1+As+v] + h h woMPtf) ' 

(85) 

which adds the missing NLO correction of the real emission 
type and also includes A$+v, representing the remaining NLO 
virtual+soft corrections. 

The important point is the definition of the variables OcJ,PJ 
and a", pf in terms of ct\ and zu in the presence of many "spec- 
tator LO gluons". An extrapolation of the one-gluon matrix 
element all over the multigluon phase space is an inevitable 
feature of any scheme combining the fixed order ME with the 
resummed ME, and there is always certain freedom in doing 
that. The above extrapolation is done in terms of zj and aj. In 
the ISR part of the sum, we proceed such that first we define 
a!- = 1 — zj and next from the evolution scale aj/ai = P'- / ot^- 
we calculate P^ (apparently we proceed as if there were no 
spectator gluons). In the FSR part we proceed similarly, i.e. 
using z J we define P" = 1 — z J ■ Next, from the evolution 
scale a 2 j/a\ — P"/ a / we calculate a". 



H. Analytical integration of DIS MC distributions and 
determining A$ +v 

A remarkable feature of the complicated multigluon dis- 
tribution defined within the exact phase space (with the full 
energy-momentum conservation) is that it can be integrated 
analytically! This integration result will help us to determine 
the NLO soft+virtual correction A$+v and will also be used in 
the numerical cross-check of the MC code. 

The result of the analytical phase space integration for the 
DIS MC reads 



2„NLO 
MC 



dtdxg 



dxidz8 XB=XlZ Di(t,xi) 



W (y B )(l +A s+v )5 l=z + C I (z,y B ) 



+ C F {z,y B )+W (y B )C [ r\z) 



(86) 



Q(zjb) = 



2C F a s 



d 3< E(k) P/(a,p,y B )9p <a 8i_ z=a 



2Cpas <d?<E{k) 



x [W(u,$,y B )-Wo(y B )P(l-a)}Q 



C F {z,y B ) 



2C F a s i 3 



d J E(t) pF(oc,p,y B )6p >a 5i_ z=a 



2C F a s 



d i T,{k) 



x [W(a ! p,y fl )-Wb(y s )P(l-P)]e p>a 8 1 _ z=a , 

(87) 



and 



cm 



2C F a s 



t / 3 £(^) J p(i-p)e a< p5 1 _3=aea> 5 



-8,=iSf(5), 
dzC l f\z)+ = 0, 

S F (5) = ^ / d^(k) P{\ - P)9 a<p 9 a>5 



1-5 



dzChz) 



(88) 



The plus prescription for (z)+ is provided by the Sudakov 
form-factor of the MC. Ci(z,y B ) and C F (z,y B ) are completely 
finite/regular, without any (...)+ parts. The IR regulator 8 will 
drop out at the end. 

Let us now find out A$+v of the MC by means of com- 



paring/matching the 1st order eq. (58 1 and/or eq. (65 i with 



eq. ( 86 1 truncated also to the 1st order. Going back for a mo- 



ment to n = 4 + 2e we find out that the 1st order bare PDF of 
the LO MC is 



£>/(/, */) I ht. ord. =8(l-xi)+K I (xi,e), [ dzKj(z,€)=Q. 

Jo 



Also, as anticipated, the contribution c3 (z) cancels the coun- 
terterm in C F (z,y B ) 



W (y B )C^(z) +C F (z,y B ) = D F {z,y B ,h) 

-W (y B )S F (8)8(l-z), 



DF(z 7 y B ,?>) 



2C F a s 



d 3 -E(k) 5i ^ocOi z>6 



(89) 



W(a,p,y B )e p; 
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Altogether the 1st order truncation of eq. ( 86 1 reads 



d 2 a NW 

U °MC(1) 

dtdxs 



\ 2 Q {5i=, B W (.y B ) [1 + A s+V - S F (8)] 
C,(z,y B )+D F (z,y Bl 8) = 2 ^ fd 3 £(k) 



x {W(a,p,y B )-Wo(^)P(l -a)e Pl<ai }5 1 _ z= „e 1 _ z>5 
[Wo(yB)B l - z >&C 2 (z)~y 2 B C L (z)l 

(90) 



where 



C s 2 (z) 



C F a s 



1 



ln- 



1 



3 1 



n I2(l-z) l-z 41-z 
Remembering that (cf. eq. §65\ ) 

cf c (z) = (c s 2 (z)) + = Q^ z>8 c s 2 (z)-8 z .. 

,■1-8 

r(5)= / dxCUx), 
Jo 

we finally get 



(91) 



° °MC(1) _ l7l( dq a QED 
dtdxg t 2 

x {8i = . VB W (y B )[l +A S+V -S F (8) + T(8)] 

+ Wo(y fl )A/(* flj e) + W (y B )C^ c (x B ) -y\ C L (x B )}. 

(92) 

Comparing with the NLO-complete (real+virtual) calculation 



(eg. eq. (63i), we see that the matching with the above MC 
implementation dictates the following relation (the Adler sum 
rule for F 2 ): 



A S+V =S F (5)-T(8) 



dz[C [ '\z)-C s 2 (z)}. 



(93) 



The above is finite in the 8 — > limit. This is not surprising, 

because o3{z) integrates the FSR counterterm, while C 2 (x) 
comes from the ISR-subtracted exact ME - they both coincide 
in the FSR collinear limit, while the ISR collinear singularity 
is already removed from C 2 (x). 

Summarizing, the complete analytical result for the struc- 
ture function from the DIS Monte Carlo (angular ordering) 
defined in eq. ([84} takes the following final form: 



a °DIS 

dtdx B 



2m 2 QED Q 2 q 



j dxdz §x B =xz Di (t , x) 



W (y B )(l+A s +v)?>i 
C F a s < 1 1 ~ 2 



Wo(y B )C? c (z)-y B C L (z 
3 



— rln(l -z) h 1 

2(l-z) V ' 41-z 



C L {z) = z. 



(94) 



The above formula is "ready to go" for numerical comparison 
with the Monte Carlo. 



The virtual+soft correction A$+v is given by eq. (93 I, more 
precisely 



As+v — 



C F a s 



W-^-It 1 - 

o t l ~ z 41-z 

.2 



1 



ln(l 



, 3 1 

■z ) + - 



3 

4^ 
-1 



C F a s 



2(1 -z) 

dz { — ln(l-z) --7--Z 



4 1-z 
1 

4~4~ 



1 



3 

2 Z 



= 0. 
(95) 



The above is just the result of the rigorous NLO calculation. 

Notice also that the MC result features in a natural way the 
exponentiation of the distributions like 



/( z ) = 5(l-z) + . 



C F a s /ln(l -z) 



C F a. s ln(l-z) c^i ln 2 (1 _ 



(l-z) 



(l-z) 



Jo 



dzf(z) = l. 



Such an exponentiation can be included in the analytical for- 
mula. 

Last but not least, let us write explicitly the difference be- 
tween the coefficient functions of the standard MS factoriza- 



tion scheme of Eq. ( |59[ > and the MC factorization scheme of 
Eq. ((94P 



AC 2 {z)=C 2 {z)-C^{z) 
C F a s f 1+z 2 



7t {2(1 -z) 



In 



(l-z) 



l-z 



(96) 



The above function should be used to correct the existing MS 
PDFs before using it to fix input in our MC. Alternatively, 



the coefficient function of Eq. ( 94 1 should be used to fit DIS 
experimental data with the PDF function compatible with the 
presented MC^ 



I. Factorization scheme independent relation between DY and 
DIS processes 

In spite of the change of the factorization scheme in 
the MC, the factorization scheme independent and the 
regularization-independent relation of AEM (eq. (91) in [30 1) 



35 Here, the usual MS assignment t = ire~' ih - is done, see refs. 13911401 . 

36 Similar corrections have to be also determined for the NLO inclusive ker- 
nels, once the NLO corrections are included in the ladder part of the MC, 
see the first incomplete results in ref. 1171. 
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should be reproduced exactly, if we claim to protect the uni- 
versality. Let us verify it. The original AEM relation reads- 3 7 | 



*f M (z)=f q ,DY-2f g , 2 



Cf oCj- 



i 



2(l-z)+ 



(-3-2z) + (l+z 



2W M 1 ^) 



1-Z 



(97) 



5 



.(-3-2z)++ (1+z 2 ) 



2(l-z) + 



1-Z 



Using the result of the analytical integration of the DIS MC, 



eq.(94i, 



Cl(z) 



1+z 2 
'2(1 -z) 



ln(l-z) 



3 1 ,3 

H + ^z 

4 1 -z 2 



and the analogous analytical result for the DY MC of eq. (44 1 



C F a s (2 2 7\ C f a, 

c 2 (z)=5 z=1 — ^ - 4 j+— [-(1- Z ) + ], 

we obtain from our two MC implementations the same result 
C 2 (z)-2C^(z)=Af M (z). (98) 



In this way we have reproduced the AEM II30II result for the 
MC factorization scheme, confirming its universality. The 
above agreement with the AEM result is easily traced back 
to the fact that it holds already for the difference of the unsub- 
tracted coefficient functiona 38 ! 



JDY W ~ 2C 2.bare W = \ W + 2 ln 



In- 
r 



(99) 



(cf. eqs. p7j ) and (58 1) and because in both MC formulae the 
same ISR counterterm of eqs. (40 » and (64i is subtracted 3 ^] 



In particular, terms due to the e^p component in the y-trace 
present in the DY and DIS coefficient functions necessarily 
cancel out. 

It is fair to mention that in ref. 11301 the relation of eq. |97]i 
is treated as the pQCD result for the coefficient function of 
the DY process in the DIS factorization scheme. On the 
other hand, this relation can be turned into an experimentally 
testable relation between the structure functions of the DY and 
DIS processes, testing the important principle of universal- 
ity (process independence) of collinear singularities in pQCD 
predictions, independently of any particular choice of the fac- 
torization scheme and the PDFs. 



37 Using again /„' ^ (1 +z 2 )( ln ^) , 



7 
4 ■ 



The last term is of course absent for the usual assignment s = t = fj 



38 

39 The ISR counterterm is defined in the DIS and DY processes at the ex- 
clusive level, involving P' and kinematic mapping, so the statement that 
"it is the same" is more non-trivial than in the case of the CFP inclusive 
counterterms. 



V. SUMMARY AND OUTLOOK 

We have presented a complete method for implementing 
NLO corrections to the hard process in the LO MC for DY 
and DIS processes. This method was originally developed for 
introducing NLO corrections in the ladder MC |[T5l[T6l . there- 
fore it is well suited to be extended to include NLO corrections 
in both hard process and ladder parts. 

The presented method is based on a new factorization 
scheme lfl6l [l7 il extending the collinear factorization theo- 
rems 1 4 26 1 to the fully exclusive (unintegrated) form, which 
can serve as a base for the MC distributions. All differences 
between the MS and this new MC scheme are kept under strict 
control, and we are elaborating on that in a quite some detail. 
In particular, we make a powerfull cross-check of the whole 
MC factorization scheme by showing (analytically) that the 
NLO MC results are reproducing the factorization scheme in- 
dependent relation of Altarelli-Ellis-Martinelli 11301 between 



the Drell-Yan and DIS processes, see eq. (98 I. 

The main practical results of this work are the multiparton 



distributions of eq. ( 34 1 and eq. ( 84 1, for the EW boson pro 



duction in hadron-hadron collision and electron-hadron deep- 
inelastic scattering, respectively, which are ready for Monte 
Carlo implementation. These distributions feature the NLO 
corrections in the hard process part and the LO pQCD evo- 
lution in two multiparton ladder parts. The NLO corrections 
to the hard process are introduced by means of a single MC 
weight on top of the LO distributions - it is, therefore, critical 
that the LO MC covers the multiparton phase space without 
any gaps or overlaps. This is achieved by means of using the 
angular ordering, which is also essential for good control of 
the soft gluon behaviour beyond LO, already in the LO MC. 
The correct soft limit also assures good behaviour of the MC 
weight; weights are positive and small (peaked near 1). For 
the weights distributions and other numerical cross-checks of 
the presented method we refer reader to ref. iPTI . 

In our opinion this work solves the main obstacles on the 
way to the NLO MC, based rigorously on the new MC factor- 
ization scheme. There are still many less important problems 
to be solved on the way to the practical level, i.e. construc- 
tion of the MC program applicable in the LHC data analysis. 
Let us signal some of these problems and their solutions, (i) 
For simplicity in our formulas we have omitted the initial PDF 
of the quark in hadron at low factorization scale Q ~ 1 GeV. 
This can be easily included in the MC. (ii) If we are aiming in 
a fully NLO MC the ladder parts have to be upgraded to the 
NLO level, and this work is already well advanced |[T4"ifT6l . 
(iii) In the presented MC scheme the QCD coupling was con- 
stant, non-running. It is quite trivial to make it running within 
the LO MC. It will be less trivial, but also profitable, to disen- 
tangle the running-coupling effect from the NLO corrections 
in the MC implementation of the NLO ladder. This problem is 
under study and will be treated in a separate publication, (iv) 
All the MC distributions presented in this work are defined for 
quarks and gluons, hence in the practical level MC code they 
will be subject to a hadronization procedure, using one of the 
existing MC tools, such as HERWIG E) or PYTHIA Q. 

Obviously, the proposed scheme of implementing the NLO 



26 



corrections to the hard process combined with the MC parton 
showers (ladders) is different from the existing ones. In sec- 
tion llHFI we comment on the differences between our scheme 
and that of the MC@NLO QD and POWHEG HI. More sys- 
tematic comparisons with these and other schemes l22l l23l 
will be done separately, at the time of the numerical MC im- 
plementation. 

The presented method of implementing the NLO correc- 
tions to the hard process does not have any principal limita- 
tions - it can be extended to more diagrams and other pro- 
cesses. However, at the practical level its application requires 
that LO parton shower provides for the full coverage (no gaps 
nor dead zones) of the hard process phase space relevant at 
NLO level. This requirement is typically not fulfilled by the 
classic MC parton showers like HERWIG or PYTHIA. It is 
not excluded that modernized version of these MCs will pro- 
vide for better phase space coverage, notably using tools de- 
veloped for the MC@NLO and POWHEG implementations. 
Otherwise, LO parton shower has to be reconstructed, for in- 
stance using scheme proposed in the present work. (This may 
turn out to be mandatory for implementing NLO corrections 
in the ladder parts of the parton shower MC.) 

Summarizing, this work represents an important step into 
a new area in the pQCD calculations for hadron colliders in 
the MC form, in which the NLO corrections are implemented 
both in the hard process and the ladder parts in a completely 
exclusive (unintegrated) way, in full compatibility with the re- 
defined, fully exclusive pQCD factorization. 



tion of the hard process 

q(P0F) + q{poii) -> r{qi) + l + {qi)+g{k), 



(A.l) 



which is the classic EW vector boson production process in 
the annihilation of the quark-antiquark pair (Drell-Yan pro- 
cess in case of y*) decaying into a lepton pair. Note that in the 
definition below we omit the distribution of the quark (anti- 
quark) in the proton. This can be always added easily in the 
MC. 

The following kinematical variables 
x= = — = 1 - a - p. 



s P 2 

P = P0F+P0B, Q = P = POF + PQB-k, 

s = P 2 , s = P 2 = Q 2 = {P-kf = P 2 -2kP, 



(A.2) 



a = 



2k ■ pob 

p2 ' 



P = 



2k ■ pop 
~~P 2 ~ 



are used in this work. The most important of them is the in- 
variant mass squared S of the produced colourless boson. 

For the emitted gluons we are using dimensionless eikonal 
phase space parametrized in terms of various variables: 



d 3 <E(k) = 



d 5 k 1 
2fc° k 2 



1 dk+ d 2 k 



7t d§ da da 2 
2 2% a a 2 
d§ d$ 



% d<i> da dQ d<i>da „ 

2 2n a P 2ti a 1 27C (3 " 



(A.3) 
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where k = (k 1 ,k 2 ) is a transverse Cartesian 2-vector (kf 
|k| 2 = sap), the Sudakov (lightcone) variables are 



k°±k\ 



2tf 



a ■■ 



P 



2k- 



Moreover, we introduce variable a = k/oc, and the conven- 
tional rapidity variable r| is defined as 



la 1, ^ 
In — = In — 
2 P 2 k- 



ln 



Multiparticle phase space is defined as 

dx n (p; PUP2 ,...,p n )=^(p-£p,) n^. 

v ,ti ' i=\ 2 p/ 



(A.4) 



Appendix A: Kinematics of EW boson production process 

Let us consider the case of a single real (not necessarily 
collinear) gluon emission, relevant to the NLO level descrip- 



The two-dimensional phase space for massless particles is 
then 



dT 2 (Q;qi,q2) = ^d£l. 



(A.5) 
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Appendix B: 1-real gluon NLO correction, analytical integration 



We are going to integrate analytically 1-real gluon NLO correction as defined in eq. (17i. The contribution from the F 
hemisphere is easily calculable: 



C 2r (x) = 



CfOL s 
K 

C F a s 



daj c/p 

o o 

(l-*)/2 



H a+p<l fc, p<a°x=l-a-|3 P <a JC=1_a 



(l-*)/2 



(1-*) 



CfOc s #(jc) 

71 (1-Jt) 



(l-*)/2 



Jp 

1-x-p 



Cpcc s 1 —x C/?a s 1 —x 



(B.l) 



(I-jO/2 (1-*) 

dp /■ jp 



p V p 



7C 2 



7C 2 



Appendix C: Inclusive NLO factorization formula for DY MC 



We are going to prove the formula of eq. ( 36 1, representing the MC with two LO ladders and the NLO-corrected hard process, 
by means of reorganizing the phase space integration of eq. ( |34| >. Let us consider the part of the total cross section of eq. p4| 
proportional to the term j e F in the MC weight of eq. (35 1. The summation and integration over the "spectator" LO gluons 
in the B part of the phase space can be easily folded into the LO PDF. What remains to be considered is the following sum of 
integrals: 



oT = / dx F dx B £ e s F j (ft^E&JiW, ^?/W) I h ( M^P^EA Gb ^ Xb) 5xF=Ui 

J bi=1 • / S<t)n 1 \i=l 71 / i€F ryZFj) 



2C f a s 

2 ■±,( i /c,;o T1 .< Tlj _ 1 

"1 = 1 7=1 17 \i=l,t£] 

d 3 T,(kj) Buj+K^jKm^fii&PFiPBWj^Fj) G B (Z,x B ) 8 Xp =YLm- 

(CI) 

The essential step in transforming each j-th term is relabelling the gluons i — > i' such that i' = i for i = 1 , 2 , . . . , j — 1 and i' = i—\ 
for z = j+ 1, hence /' = 1,2, — 1 without any gap, and finally z = y is relabelled as / = 0. Using the symmetry of the 
integrand, integrals over zc,v can be pulled out, and the sum over adjacent integration ranges of kf = ko is factorized off: 

m-1 



of° = I dx F dx B h x=XFXB £ e~ s * j ( J] ^ 3 £fe)e,,<n,, 1 2 % C ^fe') 

7 n 1= l 7 V ,'=1 K 

n\— 1 » /?] — 1 

/ d 3 T,(k ) 9i 1 ,,<T| <r|,.,_ 1 Pi («,/3f, p B \a Q ,z F o) G B (Z,x B ) 5(x f ~Zfo f~[ Zj?;' 



(C2) 



x 



i'=l J i'=\ 
The sum over the adjacent integration intervals is combined into a single integral 

Cl\ «3 a n^-\ a n\-\ 



J $idao + J $idao + J fiidao h ^ pic/ao = J piafao 



ai «2 a,i[-2 

and factorized off, while the remaining integrals over the spectator gluons i' = 1,2, ...,ti\ — 1 give rise to the LO PDF: 



n\ — 1 



of LO = Jdx F dx' F dx B ^ £ J ^n^^)^^-!^^^))^^^ 

x d 3 T,(ko)^i(S : p Fl p B \a ,z F o)G B {Z,x B )d XF=ZFOxlF 8 x=XFXB ^ C3 ^ 



■ J dx B dx' F dz F o G F {Z,x' F ) G B {Z,x B ) -C2r{ZFo) Ob(sx)8 x=XbXfZf01 
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where we have replaced the integration variable xf with zfo 
eq. ( fT7) i and evaluated in Appendix [B] eq. ( |B.l| i 



■ Xf/x' f . In the last step we were able to use the integral defined in 



The other part of the total cross section of eq. ( 34 1 proportional to the term j € B in the MC weight of eq. ( 35 i gives the same 
result. For the LO part times ( 1 + hs+v) we use eq. (31 1. 

As already noted, the key part of the above algebra is reminiscent of that in ref. [36|, except that here the resummed singular- 



ities are in the angle while in ref. [36] they are in the energy variable. 
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